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ABSTRACT 

We present a measurement of the Type I quasar luminosity function at z = 5 using a large sample 
of spectroscopically confirmed quasars selected from optical imaging data. We measure the bright 
end (Mi45o < —26) with Sloan Digital Sky Survey (SDSS) data covering 6000 deg^, then extend to 
lower luminosities (-M1450 < —24) with newly discovered, faint z ^ h quasars selected from 235 deg^ 
of deep, coadded imaging in the SDSS Stripe 82 region (the celestial equator in the Southern Galactic 
Cap). The faint sample includes 14 quasars with spectra obtained as ancillary science targets in the 
SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), and 59 quasars observed at the MMT 
and Magellan telescopes. We construct a well-defined sample of 4.7 < z < 5.1 quasars that is highly 
complete, with 73 spectroscopic identifications out of 92 candidates. Our color selection method is 
also highly efficient: of the 73 spectra obtained, 71 are high redshift quasars. These observations 
reach below the break in the luminosity function (Mj'45Q k, —26.8), although the faint-end slope is 
poorly constrained. The bright end slope is steep (/? < —3.5), with a constraint of /3 < —2.5 at 95% 
confidence. The break luminosity appears to evolve strongly at high redshift, providing an explanation 
for the flattening of the bright end slope reported previously. We find a factor of ^ 2 greater decrease 
in the number density of luminous quasars (M1450 < —26) from z = 5 to z = 6 than from z = 4 
to z = 5, suggesting a more rapid decline in quasar activity at high redshift than found in previous 
surveys. Our model for the quasar luminosity function predicts that quasars generate ~ 20% of the 
ionizing photons required to keep the universe ionized at z = 5. 
Subject headings: quasars: general 
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1. INTRODUCTION 

The number density of quasars evolves strongly with 
redshift, a conclusion reached shortly after the ini- 
tial identification of cosmological redshifts for quasars. 
Quasars i ncrea se in number with increasing redshift 
(jSchmidtl 119681 ) until z ~ 2.5, when quasar activity 
peaks, as surveys for higher redshift quasars show a 
steep decline in numb er fOsmer' 1982; Warren et al."1994'; 
Schmidt et all 119951 [Fan e t al. ,2001b; Richards et al. 



2006[) . Quasars are a ssociated with a ccretion onto super- 



massive black holes ()Salpeteiill964[) ; how the first such 
black holes grew from initial seeds and were triggered 
by gas accretion to become luminous quasars remain 
key questions about the evolut ion of the early universe 
(see, e.g., the recent review by iVolontcri 20ic | . Indeed , 
quasars are now observed to z ^ 7 (pvTortlock et al.ll2012D . 
indicating that the mechanisms that drive them were in 
place within 0.8 Gyr after the Big Bang. 

The quasar luminosity function (QLF) is one of the 
most fundamental observational probes of the growth of 
supermassive black holes over cosmic time. The QLF is 
generally found to have the form of a broken p ower law 
(|Bovle et al.l Il988l; |PeillT995l: iBovle et al.l 12000) . with a 
steep slope towards high luminosities and a flatter slope 
extending to low luminosities. Qbserved evolution in the 
shape of the QLF with redshift — e.g., a change in the 
power law slopes or the location of the break luminos- 
ity — may provide insight into the physics of black hole 
growth. Many studies have found e vidence for evolution 
of the QLF in la rge quasar s amples ([Schmidt et al."1995'; 
Fan et all 12001 b; Richard s et al.l [f006; Hopkins ct al. 



2007ht IC room et al.ll2009ll . One key consequence of this 



evolution is the observed "downsizing" of quasar activity: 
where by the spatial den sity of more luminous objects 
peak s at higher redshifts (iCowie et al.l 120031; lUeda et al.l 



[2551 IHasinger et al.ll2005l: ICroom et al.ll2009l)~ 

Theoretical models for quasars have explored a vari- 
ety of physical processes that may drive the triggering of 
quasar activity. These models make testable predictions 
about the evolution of the QLF based on the evolution of 
the triggering mechanisms. For example, models where 
quasar activity is instigated by mergers of gas-rich galax- 
ies (e . g., [H crnquist 1989; Carlbcrg 1990; Cattanc o et al.l 
1999; Kau ffmann fc Haehneltll2000t iHopkins et al.l 120061 : 
Shen ,20090 tie the QLF to the evolution of the merger 



rate of dark matter halos in cosmological simulations. Al- 
ternatively, some hydrodynamical simulations point to 
rapid inflows of cold gas along filamentary structures 
in the cosmic web a s the primary fueling mechanism 
at high redshift (e.g.. iDi Matteo et al.l[2012t ). It is also 
possible to predict the QLF wh ile being agnostic to the 
specific triggering mechanisms; iConrov fc White! (|2012[ 1 
describe a model for populating galaxies with accreting 
black holes that directly relates the evolution of quasars 
to that of their host galaxies. 

Furthermore, feedback from black hole accretion is ex- 
pected to play an important role in regulating the growth 
of black holes and the duration of quasar activity (e.g., 
IDi Matteo eral[[2005t [Hopkins et al.|[2006[ ). As the QLF 
is a convolution of the black hole mass function and the 
Eddington ratij^ distribution, these processes will alter 
the shape of the QLF (Hopkins ct al. 2005). While the 

The ratio of the mass accretion rate onto the black hole to 



low-redshift QLF is reasonably well measured by surveys 
at optical. X-ray, and mid-infrared wave lengths (see, e.g., 
the compilation of lHopkins et al.ll2007b[ ) , observations at 
high redshift are less constraining on these models, even 
though it is at high redshift where, for example, differ- 
ent feedback models make signifi cantly different predic - 
tions for the quasar population ([Hopkins et al.|[2007aD . 
Feedback-regulated models for quasar activity tend to 
predict str ong evolution of the bright-end slope with red- 
shift (e.g., fKauffmann fc Haehneltl l2000': 'Wvit he fc Loebl 
|2003| ). and also provide a physical ex planati on for the 
observed downsizing trend [Scannapieco et al.l ([2005| ). 

Measurements of the high redshift QLF also allow es- 
timates of the contribution of quasars to the ioniza- 
tion state of the intergalactic medium (IGM) during 
and after reionization. Quasars are unlikely to pro- 
duce enough ionizi ng photons to be the primary driver 
of reionization (Fan et al.ll2001al ) although current con- 
straints on the quasar ionizing photon budget are lim- 
ited to extrapolations from the brigh t end of the QLF 
and a handful of fa int z ~ 6 quasars ([Jiang et al.ll2009l : 
IWillott et alllMol ). as well as upper limits on moderate 
luminosity AGN from X-ray surveys (Barger et al. 200^; 
iFontanot et al.l[2007| ). Quasars do have a much harder 
spectrum than stellar sources, but constraints from the 
soft X-ray background limit the contribut ion of high en- 
ergy photons to large -scale reionization ([Diikstra et al.l 
120041 iMcQuinnI [20121 ). Nonetheless, there are few ob- 
servational constraints on the faint quasar population 
during the epoch of reionization. And while quasars 
may not be directly responsible for hydrogen reioniza- 
tion, they are expected to provide the high-energy pho- 
tons respon sible for He II reionization at lower red- 
shifts (e.g., Haiman fc Loebl 119981 : iMadau et all 119991 : 
[Miralda-Escude ct al. 2000). Models of IGM evolution 
thus benefit from improved observat ional constraints o n 
quasar activity at high redshift (e.g.. lBolton et al.ll2009D . 

Currently, only optical and near-IR imaging surveys 
provide the requisite area and depth to construct large 
samples of quas ars at high redshi fts (see, e.g., the com- 
pilation given in lRoss et al.l 120121 ). Throughout most of 
the 1990s, the redshift r ecord was held by a 1 9th magni- 
tude quasar at z = 4.9 ([Schneider et al.lll99lD. The fi rst 
quasar with z > 5 was discovered bv lFan "eTall (119991) m 
commis sioning data from the Sloan Digital Sky Survey 
(SDSS; lYork et al.l [20001 ). The same data were used to 
build a sa mple of 39 quasars with 3.6 < z < 5.0 to a limit 
of i = 20 ([Fan et al.ll2001bl ) and to study the evolution 
of quasars at high redshift, confirming the steep decline 
in number density at z > 3 (roughly a factor of three per 
unit redshift). The early high- redshift quasar surveys 
also found evidence for fiattcning of the bright-end slope 
relativ e to lower redshifts ([Schmidt et al .111 9 9 51 : IFan et aD 
[200iy), a result seemingly confirmed by a large, homoge- 
neous quasar sample from the SDSS extending to z = 5 
([Richards et al.l [20061 ). This type of evolutionary trend 
in the bright-end slope would contribute to a downsiz- 
ing effect, by slowing the decline in number density with 
redshift for quasars above the break luminosity. 

The SDS S-III Baryon Oscillation Spectroscopic Sur- 
vey (BOSS; lEisenstein et al.l 120111: [Dawson et al]|2012D 
aims to collect spectra of over 150,000 quasars with red- 

the maximal rate allowed by the Eddington limit. 
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shifts between 2.2 and 3.5. iRoss et all ()2012l ) present the 
QLF measured from 22,000 color-selected quasars from 
BOSS DR9, as well as two samples of variability-s elected 
quasars at z < 3.5 from ^ Palanguc-Dclabrouill e et al.l 
(|201l!) and iPalanque-Delabrouille et al.l (|2012l) . The 
BOSS spectroscopic target selection extends ~ 1.5 mag 
fainter than the SDSS, reaching below the break in the lu- 
minosity function to z ~ 3.5. Combined with the 2SLAQ 
(jCroom et al.ll200 9[), these large quasar surveys form a 
fairly complete picture of the optically unobscured (Type 
I) quasar population at z < 3. 

At higher redshifts, over 50 quasars are now known at 
z > 6. While constraints at these redshifts are weaker 
(and await future wide-area surveys for greater num- 
bers), the best determin ations to dat e indica te that the 
brigh t-end slope is steep piang et al.l 2008: W illott et al.l 
l2010f l. roughly agreeing with the z < 3.5 measure- 
ments from B OSS and the z 2.5 determination from 
iCroom et all £2009). One of the aims of this work is to 
examine the QLF at intermediate redshifts and test pre- 
vious claims for a flattening of the bright-end slope at 
z > 3. 

We present a measurement of the QLF at z — 5, com- 
bining bright quasars from the SDSS with faint quasars 
reaching nearly 2 mag deeper. These quasars are drawn 
from optical imaging data (probing the rest-frame ultra- 
violet) and thus have low intrinsic extinction, and are 
confirmed by spectroscopy to be broad emission line. 
Type I quasars. The faint sample is derived from coad- 
ded optical imaging in the SDSS Stripe 82 region. This 
imaging covers 235 deg^ to a depth more than 2 mag 
fainter than the SDSS main survey. The combination of 
medium depth and medium sky area - relative to SDSS 
and to small-area deep fields - is ideal for searching for 
rare, faint high-redshift quasars. With Stripe 82 we are 
able to reach sufficient depth to probe the faint end of 
the luminosity function, while attaining enough dynamic 
range to constrain its overall shape. This uniform sam- 
ple of z ~ 5 quasars — selected with simple criteria over 
a homogeneous imaging area — is larger than all z > 6 
surveys combined, and provides a key link in our under- 
standing of quasar evolution at high redshift. 

We first describe the coadded optical images on Stripe 
82 that provide the primary catalogs from which we se- 
lect candidates, as well as infrared imaging from UKIDSS 
and our own observations used to reduce stellar contam- 
ination (§[!]). Section 13] outlines our selection criteria, 
derived from models of quasar colors, and based on sim- 
ple color cuts. We combine spectroscopy from the SDSS 
and BOSS surveys with our own observations using the 
MMT and Magellan telescopes to build a sample of over 
70 z ~ 5 quasars to a limit of iab = 22; the observations 
and data processing are presented in section 2] and the 
catalog of quasars is provided in section [5l In section [HI 
we use our quasar color model to quantify the complete- 
ness of the survey; with this in hand we calculate the 
binned luminosity function and derive a parametric form 
for the luminosity function using the maximum likeli- 
hood technique. We explore the evolution of the QLF at 
high redshift and determine the contribution of quasars 
to the ionizing background at z 5 based on our QLF 
model. We present conclusions in section [T) 

We incorporate photometric data from several sources 



in this work; for consist ency, all magnitud es are con- 
verted to the AB system ()Oke fc Gunnlll983l ) unless oth- 
erwise note d. SDSS reports u qriz photometry on the 
asinh scale (jLupton et al.l [19991 ) . which is nearly identi- 
cal to AB at bright magnitudes. Although we use SDSS 
photometry only for bright quasars, we will note any in- 
stances where the asinh system is used. We refer to the 
SDSS DR7 imaging as "SDSS main", in contrast to the 
deep imaging we use to build our primary faint quasar 
sample, referred to as "Stripe 82" or "the coadd imaging" 
( §2.1|) . UKIDSS magnitudes are Vega-based, but we con- 
ve rt all UKIDSS r uagnit udes to AB using the values given 
in iHewett et al.l (|2006D . All magnitudes a re corrected 
for Galactic extinction (jSchlegel et al.|[l998l ) unless oth- 
erwise noted. We use a ACDM cosmology with parame- 
ters = 0.728, Vt ra = 0.272, r^b = 0.0 456, and Hq = 
70 km s-i Mpc"^ (|Komatsu et al.ll2009f ). 

2. IMAGING DATA 

2.1. Stripe 82 coadded imaging 

Over the ten year duration of the survey, the SDSS 
I/II repeatedly imaged a 2?5 x 100° stripe cen- 
tered at zero declinati on in the Southern Galactic Cap 
(|Abazaiian et al.l [20091 ). This repeat imaging increased 
in frequency during the latter part of operations as part 
of the SDSS Supernova Survey (jFrieman et al.ll2007D . By 
the end of operations, over 100 repeat imaging scans were 
obtained along Stripe 82 (its des ignation according to 
the survey geometry, iStoughtoiTet al., 2002i ). As with 
the main surv ey, the imag i ng wa s obtained with a drift- 
scan camera (iGunn et al.l |1998[) m ounted on the 2.5m 
Sloan telescope (jCunn et al.ll2006() . with nearly simulta- 
neous 54.1 s exposures in five broad optical bands (ugriz; 
Fukugita et al. 1996,). The quality restrictions applied to 
the main s urvey imaging (photometric and go od-seeing 
conditions; llvezic et al.1 120041: iHogg et al.ll200lD were re- 
laxed for the Supernova Survey in order to increase the 
temporal coverage, resulting in a greater range of data 
q uality. 

iJiang et all ()2009l ) describe deep riz images of Stripe 
82 constructed from the coaddition of 50-60 individ- 
ual imaging scans at each position on the sky. Images 
with seeing poorer than 2" or an r-band sky background 
brighter than 19.5 mag arcsec"^ were rejected. The re- 
maining images were weighted by their transparency, see- 
ing, and sky back ground noise, and coadded using the 
SWARF software (jBertin et al.ll2002t ). The final coadded 
images reach a depth roughly 2 mag fainter than single- 
epoch SDSS images, and wer e used to discove r z ^ 6 
quasars as faint as zab = 22.2 (| Jiang et al.ll2009D . These 
dat a are not the only c oadded images available for Stripe 
82: lAnnis et al.l (|2011[ ) produced coadded images made 
using a similar process but with stricter image quality 
criteria resultin g fewer ini a ging s cans used at each sky 
position, while iHuff et all ()20lH ) created image stacks 
optimized for weak lensing shear studies. 

Our start i ng po int is the riz coadded images from 
IJiang et al.l (|2009| ). which we now extend to include the 
u and g bands (Jiang et al., in preparation). The deep 
imaging in the bluer SDSS bands is being used to dis- 
cover ultra- luminous Lyman break galaxies at z > 2.5 
(Bi an et al., in preparation). In this work, we extend 
the IJiang et al.l (,2009,) search for faint, z ~ 6 quasars to 
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Fig. 1. — Median photometric uncertainty for stars with i = 22 along Stripe 82. The depth is very uniform, and with a typical uncertainty 
of o-(i) ft! 0.03 at the flux limit for our quasar survey the completeness is high. We search for quasar candidates within the region bounded 
by the white lines. The cutoff at R.A. = 20'' 32™ is due to the high stellar density along the western edge of Stripe 82, as well as the 
lack of UKIDSS coverage there. Note that the photometric uncertainties increase outside the boundary, likely due to crowding in the high 
stellar density regions. The cutoff at 5 = 0.85° is due to problems in the g-band coadded images in the northernmost two camera columns 
of the Stripe. 

lows that of the parent SDSS imaginsP^. We photomet- 
rically calibrated individual fields in the coadded cat- 
alogs by deriving zero poi nts from the ii bercalibrated 
(iPadmanabhan et al.l I2008D SDSS DR8 (|Aihara et al.l 
120 lit) photometry (derived from the best individual 
imaging scan at each location in Stripe 82), using stars 
with 15.5 < r < 20.5. We thus tie the aperture photom- 
etry from Sextractor on a field- by- field basis to PSF pho- 
tometry from SDSS. This process results in zero points 
that account for seeing and photometric variations be- 
tween the fields, but not within the fields. However, the 
large number of images contributing to the coadds tends 
to average over effects such as varying sky backgrounds 
and PSF shapes. In general, we find the photometry 
is highly consistent between fields and agrees quite well 
with SDSS photometry for the brighter objects. Figure[T] 
shows that the depth across the Stripe is highly uniform, 
and Figure[2]shows the depths reached in the ugriz bands 
in the coadded imaging. 

Quasars at z ~ 5 have similar fiuxes in the i and z 
bands. We adopt the i-band as our detection band as it 
provides greater a greater signal-to-noise ratio. However, 
at z ^ 5 the Lya emission line and the Lya forest are 
within the i-band. The fluxes in the band that defines 
our detection threshold are thus subject to the substan- 
tial variances in Lya equivalent widths, the incidence 
of strong Lya absorption systems, and variations in the 
mean IGM opacity. The z-band is less affected by these 
issues (although it does contain C IV in emission); how- 
ever, it is significantly noisier than the i-band because of 
lower CCD sensitivity and a higher sky background. We 
will discuss these issues further in Sj5l 

We further note that our survey is restricted to less 
than the full area on Stripe 82 for two reasons. At the 
time that we performed target selection, the coadds cor- 
responding to the two uppermost camera columns had 
problems with the sky background in the gi-band, ren- 
dering them inadequate for our purposes. We thus im- 
posed a cutoff of (5 < 0.85° when selecting candidates. 
In addition, the western edge of Stripe 82 (near 20*^) 
approaches the Galactic plane. Attempting color se- 

The SDSS camera consists of six col umns of CCDs th at are 
continuously exposed in a drift-scan mode l IGunn et al.ll998l) . The 
gaps between the CCD columns are filled by combining two Strips 
offset in declination (each consisting of a single scan) to form a 
Stripe. Thus the declination axis of Stripe 82 is split into 12 
columns during the imaging scans. Each scan is 13' wide in decli- 
nation, broken into 10' segments along right ascension; these image 
segments are referred to as fields. 



Fig. 2. — Stripe 82 coadded imaging depth in the five SDSS 
bands. Lines represent the median S/N measured in 1.8" aper- 
tures for isolated stellar objects (no extinction correction has been 
applied). Filled circles are drawn at the 3cr limits. The vertical 
dashed lines show the approximate r (red) and g (green) magni- 
tudes for a typical z ~ 5 quasar at our survey limit of j = 22, based 
on colors of r — j = 1.5 and g — r > 2.2. By comparison, the 95% 
completeness limits for single-epoch SDSS im aging are g 22.2, 
r « 22.2, and i ^ 21.3 l lStoughton et al.|[200^ . The depth of the 
coadded imaging in the g band allows for reliable selection of very 
red objects. 

lower redshifts (z 5), as the deep g-hand photometry 
combined with the redder bands allows reliable color se- 
lection of quasars in this redsh ift range. The dept h of our 
coadds is similar to that of the lAnnis et al.l (|2011[ ) coadds 
(see their Fig. 7); however, as described in §4.21 we en- 
countered much greater contaminati on when using the 
catalogs from the lAnnis et al.l ()2011f ) coadds for quasar 
selection. 

We use Sextractor (jBertin fc Arnout^l 1996( 1 to produce 
object catalogs from the deep ugriz images. Object de- 
tection was performed in the z-band, and fluxes and other 
measurement parameters in the other bands were derived 
from apertures centered on the i-band object position 
using the dual-image mode of Sextractor. All fluxes and 
magnitudes from the coadded Stripe 82 imaging quoted 
in this work are taken from aperture photometry using 
a 1.8" diameter. The typical seeing in the coadded im- 
ageSis 1.4", 1.3", 1.1", 1.0", and 1.1" in u, g, r, i, and 
z, respectively. 

The geometry of the Stripe 82 coadded imaging fol- 



The weightin 
better seeing data; 



scheme used when coadding the images favors 
it is u) = T X (FWHM^o-^)-!^ where T is the 
transparency, FWHM is the full width at half-maximum of the 
point-spread function (PSF), and is the variance of the sky 
background l| Jiang et al.ll2009h . 
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Fig. 3. — Example images of a typical faint quasar at z = 4.9 selected by this survey. The panels contain the griz coadded imaging, 
followed by the UKIDSS J image, and a J image from MMT/SWIRC. This object (J234601. 5-003855) is not detected in the g-band, and 
has r = 23.64 ± 0.09, i = 21.72 ± 0.02, and 2 = 21.74 ± 0.08. It is undetected in UKIDSS; the SWIRC imaging yields J = 21.39 ± 0.13. In 
the best single-epoch imaging from SDSS DR8, this object is a 9cr detection in the j-band and < Scr in r and 2. 



lection of quasars in this region would result in over- 
whelming stellar contam ination (see, e.g., Figure 1 of 
iVanden Berk et al.|[2005t ). We thus restrict our survey 
to objects with R.A. > 20*^ 32™ (this region also lies 
outside the UKIDSS imaging area described in the next 
section). The final area of the deep imaging we used 
for z ^ 5 quasar selection is 235 deg , extending from 
20^^ 32™ to 4'\ and from -1.25° to +0.85° (Fig.p. 

2.2. UKIDSS 

The UK IRT Infrared Deep Sky Survey (UKIDSS; 
ILawrenc"e~et al. 2007) consists of multiple infrared imag- 
ing sur veys with the UK IRT Wide Field Camera (WF- 
CAM; iCasali et al.ll2007l) . including a wide-area com- 
ponent (the Large Area Survey or LAS) that covers 
- 4000 deg^ to a depth of Jab « 20.6 (5cr). The LAS 
includes Stripe 82, providing shallow infrared imaging of 
our candidates and additional leverage in discriminating 
high redshift quasars from stellar contaminants. We first 
apply a loose color cut to identify potential quasar can- 
didates from the coadded optical imaging, 

{r-i)> 1.0 && (i-z) < 0.625[(r-i)-1.0] + 0.2 , 

then queried the DRS plus release in the W FCAM Sci- 
ence Archive (WSA; iHamblv et all [20081 ) for infrared 
counterparts. 

Quasars at z '-^ 5 have {i — J)ab ^ (see section |3]), 
thus at the faint limit of our survey most of our tar- 
get objects do not have a counterpart in the UKIDSS 
catalogs. However, stars with similar optical colors to 
z ~ 5 quasars (namely, M and L dwarf stars) are redder 
in the near-IR (0.5 < (« — J)ab < 1-5, see Fig. [5]) and 
thus are relatively brighter at infrared wavelengths. We 
downloaded J-band image cutouts from the WSA for all 
of the candidates pre-selected from the optical imaging, 
then performed aperture photometry using the IRAF 
aper task at the i-band position from the SDSS coadded 
imaging, measured in 2" diameter apertures. We vali- 
dated our photometry by checking against the UKIDSS 
catalog measurements for brighter objects (utilizing th e 
UKIDSS calibration as described in lHodgkin et al.ll2009[ ). 
The aperture photometry reaches J « 21.3 at 3a, suffi- 
cient to discriminate the typical stellar contaminants at 
our limit of i = 22. Thus the IR information is highly 
useful as stellar veto: of the ^ 10, 000 pre-selected candi- 
dates from the optical imaging, only ^ 650 are rejected 
by the i — J color cut we adopt ( ij3.2p using the UKIDSS 
catalog photometry, while ^ 4500 are rejected based on 
the aperture fluxes. 

2.3. MMT SWIRC 



At the time the initial candidate selection was per- 
formed, the publicly available UKIDSS data on Stripe 
82 did not extend to R.A. < 21*^ 36™. We imaged some 
of the candidates lacking J coverage on 2011 Oct 14-15 
with the MMT Smithsonian Widefield Infrared Camera 
(SWIRC: [Brown et all [20081 ) . Each object was observed 
in the J band using a 9-point dither pattern with 30s 
exposures at each position, for a total integration time 
of 4.5 min. The typical seeing was 0.6-0.9". Images 
were dark-corrected, flat-fielded, sky-subtracted, shifted, 
and combined using standard IRAF routines. Object 
photometry was measured in 2" diameter apertures us- 
ing IRAF routines, and calibrated using observations of 
UKIRT Faint Standards. 

A total of 38 objects were observed with SWIRC. This 
includes 10 objects already spectroscopically confirmed 
as quasars but that lacked UKIDSS coverage or had low 
S/N in the J-band. These objects were observed as a 
check on our J-band color selection for quasars selected 
solely by optical colors. The remaining SWIRC targets 
were taken from a sample of objects selected as high- 
redshift quasar targets based on optical colors but that 
lacked UKIDSS coverage. These objects are mainly at 
309° < aj2ooo < 360°, the region of our survey with the 
highest stellar density. Nearly all of these sources were 
readily identified as stellar contaminants from bright de- 
tections in the SWIRC imaging, and would have greatly 
reduced the purity of our spectroscopic sample had they 
not been observed with SWIRC. 

Figure [3| displays postage stamp images of a typical 
quasar within our survey. The i-band clearly has the 
highest S/N. The non-detection in UKIDSS and the 
faint detection in the deeper SWIRC imaging are ex- 
pected for a z 5 quasar; a contaminating star would 
have been easily identified from the infrared imaging. 

3. QUASAR CANDIDATE SELECTION 

Quasars are typically selected from imaging data us- 
ing their colors, but these colors are a strong function 
of redshift , and overlap with the stellar locus at certain 
redshifts (■Fanlll999HRichards et al1[200ll[200l . Optimal 
techniques have been developed for extracting the much 
rarer quasar population from the overwhelming contam- 
ination due to stars within the Galaxy. The probability 
that an object is a quasar or star can be determined by 
considering the local density of quasars and stars in the 
multidimensional space defined by the set of fluxes avail- 
able from an imaging survey. These probabil ities can be 
used to efficiently select quas ar candidates (jBovv et al.l 
120111 : iKirkpatrick et ai]|201lD . Similar techniques have 
also been applied at high redshift, where the problem is 
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Fig. 4. — Left panel: example simulated quasar spectrum at z = 4.9 (blue) compared to a composite spectrum of DR7 quasars with 
4.7 < 2 < 5.1 (gray, shifted to 2 = 4.9). The simulated quasar spectrum was randomly drawn from the simulation grid described in ^3.11 
and selected to have M1450 = —26.5, roughly the median luminosity of 2 ~ 5 DR7 quasars. The simulated Lyo forest is calculated at a 
resolution of 0.1 A; here the spectrum is rebinned to 10 A. This line of sight includes two strong DLAs (log A'^hi ~ 21.1) at 2 = 4.428 and 
2 = 4.808. The SDSS filter curves are shown as a reference for which parts of the spectrum contribute to the broadband colors. Upper 
right panel: the mean free pat h in pr oper Mpc from the simulations, represented by blue squares. The solid green line represents the 
evolution in the mfp derived bv lSongai la & Cowic (2010), along with the ztltr range (dotted lines). Lower right panel: the effective optical 
depth (tca) from the simulations, represented as blue squares with error ranges indicating the scatter among 1000 simulated lines-of-sight. 
Measurements from Songaila (2004) are shown as green triangles, and the evolution derived by J^an et al. (2006^ is sho wn by a solid 
orange line. We converged on a forest model with mfps ~ 10% lower than the mean values from ISongaila fc Cowig (120101 ): however, the 
Teff measurements are highly consistent with the observations. 

worsened by the extreme rarity of quasars and the small 
number of bands in which the candidates are detected, 
typically at low S/N (jMortlock et al.ll2012[ ). 

Fortunately, at z ~ 5 quasar colors are well separated 
from those of stars. At 2: > 4.6 the r-band is completely 
within the Lya forest and the r — i colors of quasars shift 
markedly redward of the stellar locus, while the i — z col- 
ors sample the Lya to C IV region of the rest-frame UV 
and are bluer than the stellar locus. The g-band samples 
the Lyman limit; g — r colors are generally redder than 
those of stars due to the mean forest absorption; further- 
more, most quasars will encounter a Lyman Limit System 
(LLS) near the red edge of the f;-band, leading to satu- 
rated absorption and a ij-band non-detection. Simple 
color selection at these re dshifts was used in the SDSS 
I/II ([Richards et al.l l2002D. leading to the first quasars 
discovered at z > 5 (|Fan et al. 1999) an d a total of over 
300 z > 4.6 quasars by the DR7 release ([Schneider et al.l 

In this section, we describe a model for quasar col- 
ors based on simulated quasar spectra. We use these 
model quasar colors to motivate our simple color cuts 
that achieve a high degree of both completeness and pu- 
rity. 

3.1. Simulated quasars 

iFanI (|1999D outlines a procedure for generating simu- 
lated quasar spectra using a simple empirical model for 
their spectral properties at UV/optical wavelengths. The 
simulated spectra are integrated through survey band- 
passes to generate simulated fluxes (colors) that can be 
used to define selection criteria, and to estimate their 
completeness. This requires that the spectral model re- 
liably captures the diversity of quasar spectral features, 
and thus accurately reproduces observed quasar colors. 
The basic components of the quasar model are a bro- 
ken power-law continuum, prominent UV/optical emis- 
sion lines, pseudo-continuum from Fe complexes, and 
redshift-dependent Lya forest absorption due to inter- 
vening neutral hydrogen. We have used an updated ver- 
sion of this model (described below) to reproduce the 



colors of ~ 60,000 quasar s from the SDSS-III/BOSS in 
the range 2.2 < z < 3.5 (|Ross et al.i,2012 ). We apply 
the model at higher redshift under the assumption that 
the distr ibution of quasar SEDs does not evolve with 
redshift (|Kuhn et all 120011: lYip et al.l [20P: iJiang et al.l 
12006( 1. This approach can be compared to, e.g., cloning 
the spectra of lo wer redshift q uasars to higher redshift 
(|Chiu et al.ll2005llWillott et al.l[2005l ). which carries with 
it any selection function imprinted on the lower red- 
shift sample; or using a small set of quasar templates 
(jMortlock et al.|[2012 ). which may not capture objects 
with unusual features. 

Each quasar is assigned a power law continuum with a 
break at 1100 A. The blue slope is drawn from a nor- 
mal distribution w ith /i(a) — —1.7 and a{a) — 0.3 
(jTelfer et al.l 120021 ): the distribution for the red slope 
is \x{a) = —0.5 and a{a) = 0.3. We add to this con- 
tinuum emission lines with Gaussian profiles, where the 
Gaussian parameters (wavelength, equivalent width, and 
FWHM) are drawn from normal distributions. These 
distributions are derived from fitting composite spectra 
of BOSS quasars in luminosity bins. These distribu- 
tions recover trends in the mean and scatter of the line 
parameters as a fu nction of conti nuum luminosity, e.g., 
the Baldwin Effect (fBaldwinlll977l). and blueshifted lines 
(|Gaskelll[T98l [Richards et al.ll201lD. Finallv. we include 
Fe em ission using the template of IVestergaard" fc WilkesI 
(j2001h . scaling the template in segments to match the 
Fe emission in the composite spectra. We do not include 
a contribution from quasars with unusually weak emis- 
sion lines, which accoun t for ^ 6% of quasars a,t z > 4.2 
(|Diamond-Stanic et al.| [2009'). or from Broad Absorption 
Line (BAL) quasars, which tend to have mode rately red- 
der colors (e .g., Wevmann et al. 1991; Brothe rton et al.l 
[20Mt iReichard et al.il2003D . 

The quasar model for the BOSS employs a pre- 
scription for the Ly a fores t based on the work of 
iWorseck fc Prochaskal ()201lD . This forest model is 
calibrated to observations at z < 4.6. We extend 
the model to higher redshifts by using the observed 
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Fig. 5. — Selection criteria for 2 ~ 5 quasars. The background image shows the density of isolated point-like objects on Stripe 82 with 
i < 22 and cr{u) < 2.5 (i.e., u-band dropouts); a logarithmic scaling has been applied. These objects represent the red end of the stellar 
locus and are our primary source of contamination. Overlaid are the mean color tracks for quasars from z = 4.5 to 2 = 5.5, in steps of 
Az = 0.1. The colors are derived from our simulations and include the luminosity dependence of emission line strengths; the two tracks are 
for Mi45Q = —26.7 (orange) and M1450 = —24.5 (red). The solid magenta lines are the color selection criteria used to identify 4.7 < z < 5.1 
quasars on Stripe 82. In the middle panel {riz colors), the dashed line shows the SDSS criteria used to target z > 4.6 quasars and adopted 
here for the analysis of DR7 quasars, and the dotted line shows our pre-selection criteria used to reduce the list of point-like objects for 
which we obtain aperture photometry in the UKIDSS J band. Finally, the purple points show spectroscopically confirmed 2 > 4.5 quasars 
on Stripe 82, with downward pointing triangles denoting 2 < 4.7 quasars, upward pointing triangles 2 > 5.1 quasars, and crosses quasars 
within our primary redshift range of 4.7 < 2 < 5.1. 



number densities of high column density systems from 
iSongaila Sz Cowi^ ()201Cl[ ). For the lower column density 
systems (IorNht < 17.2), we begin with the param- 
eters of iWorseck fc Prochaskal (|2011l ). and follow their 
method of deriving the number densities by simulating 
a large number of sight lines and matching the mean 
free paths (m fps) of the simulated sightlines to the ob- 
servations of Songaila fc Cowic (2010). We also force 
continuity between the column density distribution func- 
tion [dn / dNYiidz) at low and high column densities, and 
check the derived effective optica l dep th (rRff) against th e 
measurements ofl Songailal ()2004|) and lFan et all ()2006aD . 
There remains significant uncertainty in forest param- 
eters at high redshift, but our simple color selection 
method is relatively insensitive to the model for the Lya 
forest. To account for the degree of uncertainty in the 
forest absorption, we also include a "low" and "high" 
forest model, scaling the number densities of forest ab- 
sorbers up and d own by 10%, matching t he scatter in 
measurements bv ISongaila fc Cowid (|2010[ ). After com- 
parison with observed quasar colors, we find that the 
model with a mean forest density 10% greater than the 
best-fit value from the mfp measurements provides the 
best match. Figure |4] shows an example simulated quasar 
spectrum from our models, and compares the Lya forest 
observables (mfp and Toff) obtained from our simulated 
forest spectra to observations of high redshift quasars. 

3.2. Color Selection 

Based on the results from the simulated spectra, we 
concluded that the most efficient use of telescope time 
was to focus on a fairly narrow redshift range, 4.7 < 
z < 5.1, where the colors of stars and quasars are best 
separated and the selection efficiency is high. Figure [S] 
shows the gri, riz, and riJ colors of red stars and high 
redshift quasars in our sample, as well as the selection 
criteria we adopted: 

• i < 22.0 

• S/N{u) < 2.5 



• 5 -r> 1.8 OR S/N{g) < 3.0 

• r - i > 1.2 

• i-z< 0.625((r-i)-1.0) 

• i - z < 0.55 

• i-J < ((r-i)- 1.0) + 0.56 

where all magnitudes are in the AB system, measured 
within fixed apertures (as described in and have 
been corrected for G alactic extinction using the maps of 
iSchlegel et alJ (|1998[ ). S/N measurements are obtained 
from the aperture fluxes and uncertainties given by Sex- 
tractor. 

We apply no morphological criteria in our candidate 
selection after finding that the number of potentially re- 
solved objects (Sextractor CLASS_STAR < 0.8) is smaU 
(< 10% of the total selected), indicating that contami- 
nation from compact galaxies is negligible. 

These criteria are similar but not identical to the cuts 
that define the z > 4.5 inclusion regio n in the SPSS 
quasa r selection algorithm, as given in iRichards et alJ 
(j2002l ). We extend the selection nearly 2 mag fainter 
than the SDSS main survey by utilizing the coadded 
imaging on Stripe 82. We define dropout criteria in the 
u and g bands through S/N criteria within fixed aper- 
tures rather than a magnitude cut. Although the Stripe 
82 coadded imaging is highly uniform (e.g., Fig.[T]), using 
a S/N threshold more fully utilizes the u-band depth at 
a given location, as z ~ 5 quasars are not expected to 
have any detectable flux in this band. We bracket the 
redshift range of the search to z ^ 4.7 by requiring a red 
r — i color and to z < 5.1 by requiring a blue i — z color. 
Finally, we take advantage of the available infrared data 
from UKIDSS, effectively imposing a veto on objects that 
are too red in i — J. This cut is made using the aperture 
flux ratios in the i and J bands, so that objects unde- 
tected in the UKIDSS J band — as expected for z 5 
quasars — pass the flnal cut without imposing a S/N 
threshold. 
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The i-band detection catalog contains 1.5M objects. 
We apply loose pre-selection cuts in the riz colors (see 
middle panel of Fig. [5] and ^2.2\ to reduce this list to 
^ 10, 000 potential candidates for which UKIDSS images 
are downloaded and analyzed. After applying the final 
criteria listed above, we have 92 candidates to a limit of 
lAB = 22. 

4. SPECTROSCOPIC OBSERVATIONS AND DATA 
REDUCTION 

We obtain spectroscopic identifications of our can- 
didates from multiple sources by utilizing the dense 
SDSS/BOSS spectroscopic coverage of Stripe 82 com- 
bined with our own observations. First, we identify a few 
of the brightest candidates on Stri pe 82 using spectro - 
scopic data from the SDSS I/II ( Ab azaiian et al.l [20091 ). 
More recently, the BOSS collected a large number of 
spectra on Stripe 82, including several ancillary science 
programs (jPawson et al.l 120121) targeting high redshift 
quasars. Finally, we have obtained spectra for additional 
candidates (extending much fainter than the SDSS and 
BOSS spectroscopic observations) using the MMT and 
Magellan telescopes. In total, we have 73 spectroscopic 
identifications for our 92 candidates (79%). Table[l]sum- 
marizes the status of spectroscopic observations of z 5 
quasar candidates on Stripe 82, including a number of 
confirmed quasars not within the uniform sample. 

4.1. SDSS DR7 

The SDSS I/II sp ectroscopic surve y concluded with 
the DR7 release (jAbazaiian et al.l [2009). Spec- 
tra were reduced wit h the standard SDSS pipeline 
(jStoughton et al.l I2002D . Quasar target selection is de- 
scribed in [Richards et al.l J2002'); high-redshift quasars 
are targeted through two methods, both to a limit of 
i — 20.2. First, outliers from the stellar locus in griz 
space are identified as high-redshift quasar candidates, 
and second, various color cuts aim to extend the quasar 
yield within s pecific redshift ranges b eyond the locus out- 
lier selection. [Schneider et al.] ([20100 present a catalog of 
over 100,000 spectroscopically confirmed quasars drawn 
from ?a 9380 deg^ of the DR7 spectroscopic footprint 
(hereafter DR7QS0). This catalog includes 191 quasars 
at 4.7 < z < 5.1 to a limit of i < 20.2, of which 184 are 
selected by our color criteria based on their SDSS flux 
measurements. 

We utilize the DR7 data in several ways. First, we 
matched our candidate list to the DR7QS0 catalog, as 
well as to the full DR7 spectroscopic database. We iden- 
tified eight of our candidates as confirmed 4.7 < z < 5.1 
quasars in DR7QS0. There were three additional quasars 
in this redshift range that did not meet our color criteria. 
None of our candidates matched to non-quasars among 
the 1.6 million objects with spectra in DR7. 

We further employ the DR7 quasars to determine the 
bright end of the luminosity function, by constructing a 
uniform sample of DR7QS0 quasars in our redshift range 
drawn from the SDSS DR7 imaging footprint. We will 
describe this sample in i jG.l.ll 

4.2. BOSS 

The BOSS quasar survey is primarily designed to find 
quasars at 2.2 < z < 3.5 for the purpose of Lya forest 



studies ([McDonald fc Eisensteiii[[2007t [Ross et al.[[2012D . 
The Extreme Deconvolution algorithm (XDQSO) used as 
the primary quasar targeting method in BO SS is tuned 
to th is redshift range and is highly efficient ([Bow et al.l 
[20 111) . However, a large section of Stripe 82 was ob served 
by B OSS in Fah 2010 as BOSS Chunk 11 (RosseFaLl 
[2012|) : during these observations quasars were targeted 
through a combination of co lor selection (using XDQSO) 
and variability selection ([Palangue-Delabrouille et al.l 
[20 lit ) . The highest redshift quasar identified by the main 
BOSS targeting in Stripe 82, including variability, is at 
z = 4.46. 

In addition, we are conducting ancillary science pro- 
grams in BOSS that take advantage of unused fibers once 
the primary survey targets (galaxies and mid-z quasars) 
have been allocated fibers (|Dawson et al.[ [20121 ). These 
programs extend the qri and riz color cuts described in 
iRichards et al.l ([20021) to fainter quasars by utilizing ar- 
eas with multiple epochs of imaging in SDSS I/II due to 
overlapping or repeated scans. On Stripe 82 these ancil- 
lary programs utilized photometry from the ~ 20-epoch 
coaddcd images available in the DR7 CAS and described 
in Annis ct al. (201 if): these ob servations are included 
in BOSS DR9 (jAhn et al.l [201^1. We wiU analyze the 
quasars identified in overlap regions outside of Stripe 82 
in future work. 

The selection criteria fo r z > 4.6 quasar s are similar to 
the riz inclusion region in [Richards et al.] (2002), but ex- 
tended to a faint limit of i = 21.5 and using the catalogs 
from the [Annis et al.l ([201 It ) coadded imaging for tar- 
get selection. On Stripe 82, this program was allocated 
283 targets, of which 221 received spectra, but only 22 
were high-z quasars. Of those, 10 were previously known 
from SDSS DR7; thus the program yielded only 12 new 
quasars. All 12 meet our selection criteria and are in- 
cluded in this study. The ancillary program targets that 
are not quasars fall mainly into two categories: 1) ob- 
jects with fluxes contaminated by nearby bright stars, 
and 2) spurious or moving objects (e.g., asteroids). We 
matched the failed targets to the coadded image cata- 
logs and found that either they did not have matches in 
our coadded imaging, or the matches did not meet our 
selection criteria. 

BOSS spectra are collected with a fiber-fed, multi- 
object spectrograph ([Smee et al.[[2012[) and reduced with 
a pipeline described in [Bolton et al.[ ([2012f ) . We visually 
examined all of the spectra for ancillary program targets 
on Stripe 82, as well as any classified by the pipeline 
as having z > 3.6. We also cross-checked our examina- 
tions again s t the DR9 Quasar Catalog as described in 
IParis et al.l (IMl) . As with SDSS, we cross-checked our 
candidate list against all spectra in BOSS, not just con- 
firmed quasars, again finding no matches to non-quasars. 
Figure [T9[ displays the BOSS spectra of z > 4.7 quasars 
on Stripe 82. 

4.3. Magellan observations 

The remaining candidates on Stripe 82 (~ 70) were 
mostly fainter than i — 21 and required spectroscopic 
observations on larger telescopes. We first observed four 
candidates at the Magellan Clay 6.5m on 2011 Jun 11- 
13 using the Magellan Echellette spectrograph (MAGE; 
iMarshall et al.,,2008. ). MAGE provides full coverage from 
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TABLE 1 

Status of high redshift quasar 
candidates on stripe 82 



Sample 


Uniform 


All 


Candidates 


92 




... spec, ids 


73 


106 


... quasars 


71 


84 


... 4.7 < 2 < 5.1 


52 


64 



3100A to 1/im at a resolution of ~ 5800 using the 0.7" 
sUt. The typical seeing during this run was 0.6". 

Magellan spe ctra were reduced using MASE 
(jBochanski et al.l [2009), an IDL-based pipeline de- 
signed for MAGE data. Wavelength calibration was 
provided by ThAr lamps observed shortly after the 
targets and at a similar airmass, and the standard star 
Feige 110 was used for flux calibration. Of the four 
targets observed at Magellan, three are z ~ 5 quasars 
(Figure [501 one object had been previously observed in 
BOSS), and one is a star. 

4.4. MMT observations 

The bulk of our spectroscopic observations occurred at 
the MMT 6.5m telescope using the Red Channel spec- 
trograph. We used the 270 mm~^ grating centered at 
7500 providing coverage from 5500A to 9700A. We 
used either the 1" or 1.5" slit based on the seeing, pro- 
viding resolutions of i? ^ 640 and R ~ 430, respectively. 

Observations were conducted on 2011 Jun 23-24, 2011 
Oct 1-4, and 2012 May 27-28. During 2011 June nine 
objects were observed in poor seeing (1.5-3"); five were 
confirmed as quasars. Conditions during the 2011 Octo- 
ber run were fair with < 1" seeing; however, most of the 
run was lost to thick clouds and high humidity, particu- 
larly during the early part of the night. As a result, we 
obtained few spectra for candidates at 20*^ < R.A. < 22*^. 
In total, 52 targets were observed in October, of which 
41 were confirmed quasars. In the 2012 May run, 17 
additional candidates with 20*^ < R.A. < 24*^ were ob- 
served, resulting in 11 new high-redshift quasars. Condi- 
tions during this run were excellent, with 0.7-1.2" seeing 
throughout. In general we valued efficiency over quality. 
Therefore the exposure times were short, between 5 and 
15 minutes with typically a single exposure per target, 
and many of the spectra have a low signal-to-noise ratio. 
However, the quasars among the observed targets are 
easily identified by their prominent emission lines, while 
the few contaminants can be ruled out as high-redshift 
quasars by the lack of either emission lines or a strong 
spectral break towards blue wavelengths. 

Data were processed using standard longslit reduction 
techniques through a combination of PyratJ and python 
routines, including bias subtraction, pixel level correc- 
tions from flat fields generated from internal lamps, and 
sky subtraction using a polynomial background fit along 
the slit direction. Cosmic ra ys were idntified and masked 
using the LACOS routines (Ivan Dokkuml[200l . Wave- 
length calibration was provided from an internal HeN- 

■^^ Some objects were observed at a redder setting during the 
course of a program targeting higher redshift objects. 

Pyraf is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 



eAr lamp, and then corrected on a per-image basis us- 
ing night sky lines ( primarily the OH line list given by 
iRousselot et al.ir2000 ) . Standard stars were observed be- 
tween one and three times per night and used for flux 
calibration. However, the conditions were highly vari- 
able and the absolute flux calibration of the spectra is 
not reliable. Figures W] and show the MMT spec- 
tra obtained for high redshift quasars on Stripe 82. 

5. QUASAR CATALOG 

Table U] consists of our full Stripe 82 z ~ 5 quasar 
catalog, consisting of: 

• 11 SDSS DR7 quasars with z > 4.7, 

• 14 BOSS DR9 quasars with z > 4.7, 

• 59 quasars with spectra obtained at MMT and 
Magellan. 

For all objects, we provide photometry in the griz 
bands derived from the coadded imaging that formed the 
basis for our target selection. We also provide J-band 
photometry obtained from our aperture photometry of 
the UKIDSS DR8plus images, or, when available, imag- 
ing from MMT SWIRC. The table includes all z > 4.7 
quasars on Stripe 82 from the three data sources, but not 
all of these quasars are included in the uniform sample 
used to calculate the QLF (fj6]). In i i6.1.2l we will derive 
the selection probability for each quasar, this value is in- 
cluded in the catalog and we flag quasars that are not 
part of the uniform sample by assigning them a value of 
-1. 

The only line widely available at a reasonable S/N 
in our spectra is the Lya line. Our spectra do gener- 
ally cover the C IV emission region; however, this line 
shows offsets from the systemic redshift that are corre- 
lated with properties such as lum inosity and radio loud- 
ness (e.g.. [Richards et all 120 lit ). Therefore, we assign 
redshifts based on a combination of fitting the Lya linep^. 
measuring the onset of Lya forest absorption, and vi- 
sually matching a quasar template spectrum to the ob- 
served spectra. In general, the redshifts have an uncer- 
tainty of Az ^ 0.02, which is sufficiently accurate for 
calculation of a luminosity function. All of our quasars 
are consistent with a Type I classification based on their 
broad line widths {FWHM > 1000 km s'^). 

We use the simulated spectra derived from our quasar 
model f §3.ip to derive fc-corrections as a function of both 
redshift and luminosity. We derive the average correction 
for the observed i-band flux to the monochromatic lumi- 
nosity at rest- frame 1450 A (M1450) for a large number 
of simulated quasars in narrow bins of (M, z) (see ij6.ip . 
then interpolate this grid to derive an individual quasar 
correction. Figure |6] shows the luminosity-dependent fc- 
corrections derived from our quasar model and used in 
this work. Accounting for the average emission line con- 
tribution to the fc-correction as a function of luminosity 
alleviates some of the issues arising from the fact that 
our best photometry is in the «-band, which contains the 

■^^ The Lya line can also have systematic offsets, in particular 
a ~ 500 km s~^ redshift due to abs orption of the blue wing by 
the Lya forest (e.g., Shen e t al.|[2007l) . We do not account for this 
offset as it is much smaller than the uncertainties we assume for 
our redshift designations. 
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Fig. 6. — Luminosity-dependent ^-correction as determined from 
the quasar spectral simulations. The lines show the fc-correction 
from observed i-band magnitude to Afi45o for M1450 = —28 (top) 
to M1450 = -24 (bottom), in steps of AM = 1. The trend 
with redshift is fairly smooth; although Lyo enters the i band 
at z 4.6, Lycf forest absorption tends to offset its effect on the 
fc-correction. The vertical error bars at the top of the plot de- 
note the representative scatter in the fc-correction at the median 
luminosity (Afi45o = —26) at three redshifts. This scatter is also 
determined from the simulations and accounts for the variations 
in emission line strengths, continuum shapes, and Lya forest ab- 
sorption included in the quasar spectral model. These corrections 
are from the Stripe 82 simulations; the DR7 simulations produce 
similar results but are adjusted for the asinh magnitude scale. 

Lya line|^ This approach also corrects for some of the 
bias introduced by luminosity-dependent line emission 
and its effect on broadband photometric data. However, 
the distribution of intrinsic quasar SEDs at a given lu- 
minosity and redshift is quite broad, introducing scatter 
into our absolute magnitude calculations (also shown in 
Figure ini). 

Alternatively, the fc-correction c an be derived directly 
from the spectral data (see, e.g., iGlikman et al.l [20T1I 
for a discussion of spectral fc-corrections for high— z 
quasars). However, we are again limited by the low 
SjN in the continuum of our spectra. In addition, we 
do not attempt an accurate flux calibration of our spec- 
tra (indeed, the BOSS quasar spectra are kn own to have 
flux calibration errors in some instances, see iParis et al.l 
I20TI) . Thus we would need to calibrate the observed 
(noisy) spectra with the broadband fluxes from the imag- 
ing. We chose to use a template fc-correction based on 
photometry as it can be more consistently applied. 

Figure [7] shows the distribution in redshift and lumi- 
nosity for both the Stripe 82 and SDSS main samples, 
after applying our fc-corrections. 

5.1. Notes on individual objects 

In this section, we note objects that have radio de- 
tections, uncertain identifications from the spectroscopy, 
and unusual spectral features. 

J221941. 90+001256.2 (z = 4.30), 

J224524.27+002414.2 (z = 5.16): These two objects 
have radio counterparts at 1 . 4 GHz from VLA imaging of 
Stripe 82 (IHodee et al. l lMl . They are the only sources 
with counterparts in that catalog, which is derived from 
imaging over 92 deg^ to a depth of 52 /zJy beam~^. 
J221941.90-f 001256.2 has a peak flux density of Fi.4,pk = 

Note that the i-band measurement is from the coadded imag- 
ing, and is thus an average of a decade of individual measurements, 
smoothing over the variable lightcurve of each quasar. 
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Fig. 7. — Distribution of high-z quasars in our sample in lumi- 
nosity and redshift. Green triangles represent DR7 quasars, blue 
circles denote Stripe 82 quasars. The grid lines show the bin edges 
used for the calculation of the binned QLF, with dashed lines for 
Stripe 82 and dotted lines for DR7. 

0.92 ± 0.07 mJy beam^i and J224524.27+002414.2 
has Fi.4,pk = 1.09 ± 0.06 mJy beam"!. Both 
sources also have counterparts in the Faint Images 
of the Radi o Sky a t Twe nty- Centimeters (FIRST) 
catalog (Bec ker et al.l |1995[ ). with peak flux den- 
sities of i^i.4,pk = 0.87 ± 0.10 mJy beam~^ and 
-F'i.4,pk = 0.92 ± 0.10 mJy beam~^, respectively. Neither 
is included in our uniform sample as they lie outside the 
deflned redshift range. J221941.904-001256.2 also shows 
strong broad absorption line (BAL) features. 

J223907.56+003022.6 (z = 5.09), J232741.35- 
002803.9 (z = 4.75), J021043. 16-001818.4 (z = 
5.05): These three objects also have FIRST counter- 
parts, with peak flux densities of Fi.4,pk — 1.35 ± 
0.10 mJy beam-i, Fi.4,pk = 1.24 ± 0.12 mJy beam"!, 
and i^i.4,pk — 6.82 ± 0.10 mJy bearn"^, respectively. All 
three are included in the uniform sample. 

J211158. 01+005302.6 (z = 4.98), 

J234730.56+002306.3 (z = 4.71): We identified 
these objects as quasars based on their discovery 
spectra; however, the spectra are noisy and the identi- 
fications were uncertain. We later confirmed them as 
quasars with MMT observations on 2012 Aug 25. 

J211225. 39-000141.3 (z = 4.67), J030315.05- 
000347.6 (z = 4.72): These spectra also have low S/N. 
Both objects appear to have a Lyman break feature 
(more evident in the 2D spectra) and Lya and N V emis- 
sion features, with apparent N V absorption troughs. We 
include these objects as confirmed quasars. J030315.05- 
000347.6 may have BAL features. 

J222018.48-010146.8 (z = 5.62), J000552.33- 
000655.6 (z = 5.86): These two quasars were targeted 
as part of an SDSS+UKIDSS high-z BOSS quasar ancil- 
lary program. J 000552. 33-000655. 6 was first reported in 
IFan et all (|200l . 

6. RESULTS 

6.1. Survey Completeness 

We use the simulations described in section |3T] to es- 
timate the completeness of our selection criteria. To de- 
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Fig. 8. — Spectroscopic completeness of z ~ 5 quasar candidates 
from the uniform sample on Stripe 82 (jab < 22). The histogram is 
divided into objects without spectroscopic observations (gray), and 
those with spectra obtained through DR7 (green), BOSS (cyan), 
and MMT/Magellan observations (blue). The orange line rep- 
resents the spectroscopic incompleteness correction applied when 
calculating the QLF in order to account for candidates without 
spectra. 

rive a selection function, we construct a grid of simulated 
quasars distributed evenly in (M1450, z) space with 200 
quasars per bin of AM = 0.1, Az = 0.05. Each object is 
sampled from models for quasar emission properties as 
outlined in § 13.11 The simulated quasars are then passed 
through the color cuts after adding photometric errors, 
and the fraction of objects within each bin that pass the 
cuts provides an estimate of the completeness for a given 
luminosity and redshift, under the assumption that our 
quasar model accurately represents the intrinsic distri- 
butions of quasar properties. 

We generate two simulation grids. The first adopts 
photometric uncertainties typical of the SDSS main sur- 
vey and uses asinh magnitudes; we use this simulation 
to derive the completeness of the SDSS quasar survey 
at z ~ 5. The second uses photometric uncertainties 
from the coadded imaging and is used for the Stripe 82 
sample. The uncertainties are determined by fitting a re- 
lation to the uncertainties of stellar objects in the Stripe 
82 catalog as a function of band flux. 

6.1.1. DR7 Completeness 

Th e SDSS quasar selection algorithm (jRichards et al.l 
[200I ) was not finalized until the survey was already in 
progress, thus to construct a st atistical sample f rom DR7 
we follow the method given in iRichards et al.l ()2G06) to 
identify DR7 quasars from regions with uniform target 
selection. This limits the final area to 6222 deg^. We 
further cut the DR7 sample to only objects selected by 
the riz color inclusion regions (excluding objects selected 
only as griz stellar locus outliers); this results in the 
loss of only a few objects, but greatly simplifies the cal- 
culation of the selection function. After restricting to 
the uniform targeting area and applying the color cuts, 
DR7QS0 provides 146 quasars with 4.7 < z < 5.0. 

We use the simulated quasar photometry combined 
with the riz color cuts of lRichards et al.l ()2002l ) to derive 
the selection function for z ^ 5 quasars in the main SDSS 
sample (Figure[9l). We add a 5% correction for photomet- 

We do not double-count DR7 quasars that lie on Stripe 82 
when calculating the LF, as the uniform targeting area does not 
include Stripe 82. 



ric incompleteness (i.e., objects lost due to crowding or 
proxi mity to bright stars or galaxies, see IRichards et all 
l2006f ) and spectroscopic incompleteness of 5% (objects 
selected by the targeting algorithm but without a spec- 
trum in DR7) . We determined the latter incompleteness 
by querying the DR7 Target database for objects with 
the QSO_HIZ fiag, finding that 5% lack spectra. The 
missing spectra are mainly due to fiber collisions that 
result when tiling spectroscopic targets, due to the re- 
striction that fibers on a single plate must be separated 
by > 55". 

6.1.2. Stripe 82 Completeness 

Although the Stripe 82 selection function is obtained 
from the same quasar model as for DR7, it is evident from 
Figure IH] that the Stripe 82 selection probes a different 
population. First, the Stripe 82 data are much deeper 
and thus the completeness remains high to lower lumi- 
nosities. Also, the color selection criteria are somewhat 
different. Figure [10] compares the selection efficiencies 
for DR7 and Stripe 82 as a function of redshift. It is 
noteworthy that the selection function value is higher for 
fainter objects at some redshifts. This is because the 
selection function calculation averages over a range of 
colors in each (M, z) bin, with the colors depending on 
the continuum slopes, line strengths, and other features 
sampled from the model. At z ~ 5, our color selection is 
very sensitive to the Lya equivalent width, as shown in 
Figure [TTJ Because the EW increases at lower luminosi- 
ties and Lya is fully within the i-band in our redshift 
range, the colors of fainter objects are redder in r — j 
and bluer in i — z than brighter objects, increasing their 
likelihood of meeting our selection criteria (see also Fig- 
ure [5]). This effect is not captured by models that do not 
account for the Baldwin Effect. 

We estimate the photometric coni pleteness to be 95% 
(see, e.g. Fig. 7 of lAnnis et al.ll2011|) . The spectroscopic 
coverage of the Stripe 82 candidates is high: ~ 90% of 
candidates with iab < 21.5 have spectra. Near the faint 
limit, the completeness is a bit lower: 69% of candidates 
with 21.5 < lAB < 21.7 and 63% of candidates with 
21.7 < ZAB < 22 have spectra. Figure H] shows our spec- 
troscopic completeness as a function of observed magni- 
tude, as well as the correction we use to account for the 
missing spectra in our luminosity function calculation. 
In total, 19 candidates do not have spectroscopic iden- 
tifications; by applying this correction we assume that 
they are a random subsample of the full candidate set. 
This assumption is fair given that the choice of which tar- 
gets to observe was largely constrained by sky location 
and weather conditions, not properties of the candidates 
themselves (such as color). 

The Stripe 82 survey includes 52 quasars at 4.7 < z < 
5.1. The Stripe 82 sample is not only highly complete but 
also highly pure: out of 73 candidates with spectroscopic 
identifications, 71 are high- redshift quasars (z > 4), and 
52 (71%) are quasars in the targeted redshift range. For 
completeness, we provide in Table [5] a list of the candi- 
dates that were either not observed spectroscopically (19 
objects), or were found not to be high- redshift quasars 
(2 objects, both have no emission features and appear to 
be stellar continua). 

6.2. Binned Luminosity Function 
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Fig. 9. — Selection functions for the SDSS main sample (left, iasinh < 20.2) and Stripe 82 (right, jab < 22). Note the difference in scale 
on the y-axis, as the Stripe 82 survey reaches nearly two magnitudes fainter than the SDSS main. The selection probability is calculated 
in each {M, z) bin by determining the fraction of the 200 simulated quasars in each bin that pass the selection criteria, averaging over the 
SED distribution. 




Fig. 10. — Comparison of the selection functions for the SDSS 
main sample (green) and Stripe 82 (blue) . The lines represent slices 
at fixed observed magnitude through the functions shown in Fig- 
ure|9] For the SDSS main, the slices are at iasinh = (19,20,20.2), 
and for Stripe 82 they are at jab = (19.1, 20.2, 21, 21.5, 21.9). The 
lines become thinner and lighter with increasing magnitude. We 
restrict our uniform sample to 4.7 < z < 5.1 where the complete- 
ness is highest (gray shaded region). Although the completeness 
is relatively high at z ~ 4.65, the selection function is steep here, 
which would make incompleteness corrections very sensitive to sys- 
tematics in the selection function model. 

We first calculate the luminosity function from the 
combined SDSS main and Stripe 82 z 5 quasar samples 
by dividing the sample into discrete bins of luminosity 
and redshift. Guided by our completeness calculations, 
we restrict the sample to the interval 4.7 < z < 5.1, 
where the selection efficiency is relatively high. We re- 
fer to this as the uniform sample. We use a single red- 
shift bin, ignoring any evolution of the QLF parameters 
over the width of th e bin. In particul ar, the redshift 
evolution derived by iFan et aTl (j2001 b^ predicts a de- 
cline by a factor of ~ 1.5 in space density from z — A.? 
to z — 5.1. We do not account for this evolution in 
the binned QLF, but it will be incorporated below in 
a maximum likelihood fit to each quasar. We calculate 
the binned luminosity function using the l/Va method 
(jSchmidt 1 968; Avni fc Bah call 1980), including the cor- 
rection of iPage fc Carreral |2000l ) . The calculation is 
performed separately on the DR7 and Stripe 82 data, 
accounting for the differences in sky area, depth, and 
selection criteria between the two samples. 

Table [2] provides the binned QLF for both the SDSS 
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Fig. 11. — Dependence of the Stripe 82 selection function value 
derived from the simulations on the Lya equivalent width. The 
simulated quasars are divided into three bins of Lya EW, with 
median values shown on the plot. The anticorrelation between 
EW and luminosity means that objects with greater Lya EW are 
fainter on average, and thus have a lower selection probability. On 
the other hand, the Lya line itself introduces redshift-dependent 
trends. For example, objects at z > 5 with strong Ly« emission will 
have redder r — i and bluer i — z colors, and thus be more likely to 
pass our selection criteria. This effect is apparent in Fig. (5] which 
shows that lower luminosity sources fall within the color selection 
criteria over a wider redshift range than higher luminosity sources. 



main (DR7) and Stripe 82 samples. We include the num- 
ber counts in each magnitude bin, as well as the corrected 
number counts after accounting for all sources of incom- 
pleteness. The binned QLF data is also displayed in Fig- 
ure [T^l The SDSS data show a steep drop in the number 
density at the bright end; from the combined data it is ev- 
ident that the QLF becomes shallower towards lower lu- 
minosities. Fitting a single power law to the binned data 
from both surveys at the bright end (M1450 < —27.0) 
r esults in a stee p slope of /3 = —3.7. 

iShen fc Kellvl ()2012f ) have also calculated the binned 
QLF of S PSS quasars at z = 4 .75, repeating the DR3 
analysis oflRichards et al.' (|2006') for the large r DR7 sam- 
ple. O ur methodology differs from that of Shen fc Kellvl 
(j2012| ) in several respects. Most notably, we have recalcu- 
lated the selection function and fc-corrections using our 
new quasar model. In addition, we restrict the quasar 
sample to color-selected objects (excludiiig those identi- 
fied only as stellar locus outliers). The [Richards et all 
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Binned QLF 
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^ <& is in units of Mpc mag '-. 
^ is in units of 10^^ Mpc~'^ mag^^. 

(|2006f) selection function has a vahie of 1 at z > 5, 
much higher than the vahies we obtain (see Fig. [TU)) . 
and inconsistent with our previous finding f ijG.l.ip that 
few quasars were selected by locus outlier criteria alone. 
Howev er, although our select ion function disagrees with 
that of [Richards et al.l H2006f ) by as much as a factor of 
^ 2 at z ^ 5, the highest r edshi ft bin (4.5 < z < 5.0) in 
both [Richards et all (l2006h and IShen fc Kellvl (IMl) is 
dominated by objects near the low redshift edge of the 
bin. Finally, we use a slightly different method for calcu- 
lating the spatial area of the survey; however, we obtain 
a similar result (6222 deg^ vs. 6248 deg^) Figu re 
compares our binned QLF to [Richards et al.l (|2006D and 
iShen fc Kellvl (|2012D ; the agreement is generally good, 
though differences of ~ 20-30% may still be attributed 
to the different approaches used. 

6.3. Parameter Estimation From Maximum Likelihood 

We now derive parametric fits to the observed data us- 
ing maximum likelihood estimation. The maximum like- 
lihood estimate for a luminosity function $(M, z) corre- 
sponds to the minimum of the log likelihood function 

5* = -2^ ln[$(M„ Zi)]+2 / / $(M, z)p{M, z)—dMdz 

i 

where the first sum is over the observed quasars, and the 
second is over the full luminosity and reds hift range of the 
samp le a nd pr ovides the normalization (Marshal l et al.l 
119831 see lFan e t al. 2001a and Kelly et al. 2008 for alter- 
native derivations of the likelihood function). The term 
p{M, z) is the probability that a quasar with absolute 
magnitude M (for our purposes, M1450) and redshift z 
is included in the survey; i.e., the selection function as 
derived in section 16.11 including all sources of incom- 
pleteness. Confidence intervals are derived from the like- 
lihood function using a distribution in AS — S ~ S'min 
(|Lampton et al.lll976[ ). 

Over a wide redshift range, the quasar luminosity 
function is found to be well fit by a double power law 
(|Bovle et al.|[T988l) , 
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where M* is the characteristic luminosity at which the 
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Fig. 12.— Binned QLF at (z) = 4.9 from SDSS main (blue) 
and Stripe 82 (red). The dotted line shows a single power law fit 
to points with Afi45o < —27.0 (/3 = —3.7). The departure from 
a single power law is evident when the faint Stripe 82 data are 
included. For comparison, we show previous calculations of the 
binned QLF at z = 4.75, corrected for ou r cosmology and s hifted 
to z = 4.9 using the redshift evolution from' Fan et al.l Il2001bl 'l. The 
orange points represent the SDSS DR3 calculation from an area of 
1622 deg2 I Richards et al.il200Ql, while the purple points show the 
SDSS DR7 calculation from IShen fc Kellvl l l20l3 'l using a sample 
nearly identical to our SDSS main sample covering ~ 6200 deg^. 
The dashed lines co rres pond to the Hopkins et al. (2007b) QLF 
model discussed in i|6.4l most of the constraint on this model at 
2 = 5 comes from the iRichards erill l(2006l ') points. 

function changes slope from steep at the bright end (/3) to 
shallow at the faint end (a). The four QLF parameters 
may evolve with redshift, possibly in an interdependent 
manner. Given the limited redshift range of our survey, 
we will only account for the steep dec line in nuinber den - 
sity at high redshift using the fit of iFan et al.l ()2001b[) : 
<^>*{z) = <^*{z = 6) X lO'^(^-s), with k = -0.47B 

Even with a sample of nearly 200 quasars spanning 
AM « 4, there are substantial degeneracies in fitting 
the data with this parameterization. In particular, while 
the departure from a single power law form towards faint 
luminosities is clear from the binned data, the faint end 
slope itself is not well constrained. This results in consid- 
erable freedom in the placement of the break luminosity. 
We thus perform several fits while fixing one or more 
parameters. 

Table[3]shows the results of several model fits with vari- 
ous choices for the fixed parameters, including uncertain- 
ties derived by varying a single parameter and calculating 
the likelihood for the best-fit solution for the other pa- 
rameters. First, we compare fits with different fixed val- 
ues for the faint end slop e. Surveys of faint z S, 3 quasars 
typically fi nd a = —1.5 (Hunt et al. 2004V Richards et al.l 
12005; Si ana et al.l [2b08; Croom ct al. 2009} ), while obser- 
vatio ns at high redshi ft favor a steeper va l ue of a = —1.7 
(iGlikman et all 120101 : llkeda et all 120111 : IMasters et all 
I2012D . albeit with large uncertainties due to the difficul- 
ties of assembling large samples of faint quasars at high 
redshift. The model fits in Table [3] are compared to the 
binned QLF in Figure [T3l 

We find that the steeper value for the faint end slope 
is marginally favored by our data, but at only the ~ 33% 
confidence level. The a = —1.7 fit has a higher break 
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Fig. 13. — The double power law fits from maximum likelihood 
analysis (Table (3)1 plotted against the binned QLF data. Each 
model fit is labeled by the choice of fixed parameters for that fit. 



TABLE 3 

MLE FIT PARAMETERS 



log4>*(2 = 6) 


-'^-'^1450 


a 


/3 


AS 




-8.ioi«:^? 

-8.46l«:i? 


-26.8lir3^ 

-^3-00-0.72 

ZD.OO_Q 3Q 


-1.70 
-1.50 
-1.70 


O ir, + 0.40 

-3.50 


3.40 
0.38 


0.33 
0.83 


-8.40lH^ 


-26.39 


-1.80 


-3.26 


0.40 


0.52 



Note. — Parameters without uncertainty ranges are fixed during 
the maximum likelihood fitting. The last two columns give the con- 
fidence level for the best-fit model (first row) compared to various al- 
ternative choices for the fixed parameters, in terms of AS = 5— S',„in, 
which follows a distribution. 

log'I>*(^) = log** (2 = 6) -I- k{z - 6), with k = -0.47. 



luminosity {M^^^q = —26.8) and steeper bright end slope 
(/3 = —4.0) compared to the a — —1.5 fit {M^^^^q = 
—25.9, (3 — —3.1). We further explore the parameter 
degeneracies by calculating the joint likelihood ranges 
for f3 and Mi^^q for the two fixed values of the faint 
end slope, as shown in Figure 1141 For a — —1.7, we 
derive 95% confidence intervals of /3 < —2.67 and 27.35 < 
-^1450 < —24.87; for a = —1.5, we derive 95% confidence 
intervals of /3 < -2.51 and -26.94 < < -23.48. 

The strong constraints on the steepness of the bright end 
slope show that a flattening o f the bright end slope at 
high redshift, as found by, e.g.. iRichards et all ()2006[) . is 
not in agree ment with our data. This is likely due to 
the fact that IRichards et al.l (|2006[) fit a single power law 
slope to the SDSS data under the assumption that the 
break luminosity was well below their flux limit; we will 
discuss this further in the following section. 

6.4. Evolution of the QLF 

We provide some context for these fits by coniparing 
to results at lower and higher redshift. iRoss et all (|2012[ ) 
recently presented a QLF measurement at 2.2 < z < 3.5 
from the BOSS DR9 quasar sample. Over this redshift 
range, the QLF is well-fit by a Luminosity Evolution and 
Density Evolution (LEDE) model, where the power law 
slopes have fixed values and the normalization and break 
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Fig. 14. — Likelihood contours for /3 and Mj^jq, for two different 
(fixed) values for the faint-end slope a. Contours are drawn at 
the 68%, 95%, and 99.6% confidence intervals, with decreasing line 
thickness. The best-fit values for each parameter are indicated with 
a cross. 

luminosity evolve in a log-linear fashion. Specifically, 

log[$*(z)]=log[<i>*(z = 2.2)]-)-ci(z-2.2) , (1) 
M*2(z) = M*2(z = 2.2) + C2(z-2.2) . (2) 

In equation [2 M^,2 = M,{z = 2) = Mu.'sn - 1-486 is 
the absolute i-band magnitude at z = 2 ([Richards et al.l 
I2006f ). corresponding to rest-frame ~ 2600 A and assum- 
ing a spectral index of = —0.5 (/^ oc i/"). The evo- 
lution in M* and <i>* implied by this model is shown in 
Figure [151 extrapolated to higher redshift to compare to 
our data. While evolution in the power law slopes is 
not well constrained by the data, our results are con- 
sistent with the values obtained from the BOSS data, 
with some indication that t he faint end slope steepens to- 
ward higher redshift [from lRoss et al.|[2012l a{2.2 < z < 
3.5) = -lA2toil and ^(2.2 < z < 3.5) = -3.53t[J;^^]. 
To further facilitate comparison to the BOSS results, we 
include a fit with the bright end slope fixed to /? = —3.5 
while retaining a steeper faint end slope of a = —1.7. 
This choice of slopes results in a fit that is fully consis- 
tent with our best fit (the indicates this model can 
only be ruled out when compared to the best-fit model 
at the 17% confidence level). However, the values for M* 
and $* do not agree with the simple extrapolation of the 
LEDE model from lower redshift. This can be clearly 
seen in Figure [TH which shows the LEDE prediction at 
z = 4.9 significantly overestimates our measurements. 
We will discuss a modified form to the LEDE model that 
provides a better fit to the high redshift evolution in H6.6[ 
At higher redshift, we compare to the z ^ 6 QLF 
measurement reported by iWillott et al.l ([20101 ) . Table [3] 
shows a fit where all the parame ters except ^* have b een 
fixed to the best-fit values from IWillott et al.l ([2010f ) for 
a fixed faint end slope of a = —1.8 (we use this value as 
it provides a better fit to our data than their a = —1.5 
fit). This model differs from our best- fit model by a 
likelihood of AS — 0.4 and is disfavored at only the 
48% confidence level. Thus the shape of our QLF at 
z 5 is consistent with the current best measurements 
at z ^ 6. However, the normalization is quite differ- 
ent: it is a factor of ~ 1 .6 higher than predicted by 
evolving the IWillott et all ((2OIQ1) model to z = 4.9 us- 
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Fig. 15. — Evolution of the QLF normalization ($*) and break 
luminosity (Mj^jq) between 2 ~ 2 and z ~ 5. The points at 2.2 < 
z < 3.5 come from the BOSS DR9 QLF IIRoss et al.|[20Tg) . using a 
sample of ~ 6000 variability-selected quasars from Stripe 82. The 
point at z = 4. 9 is from the best-fit model in Table [3] The point 
at 2 = 4 is fromlMasters et al.l II2012I ) and the one at 2 = 6 is from 
IWillott et anil20ia ). using their a = —1.8 fit (uncertainties for the 
parameters were not reported for this fit). All points have been 
corrected to match our cosmology. The log-linear LEDE model fit 
to the BOSS data is shown as a dark blue dashed line; the modified 
form of this model discussed in ^6.6l is shown as a light blue dashed 
line. 

ing k = —0.47, as they adopted for their fit (they ob- 
tained log$*(z = 6) = —8.6, in contrast to the value of 
log$*(z = 6) = —8.4 we obtain when fitting the shape 
of their QLF to our z — 5 data) . Figure [12] illustrates 
the discrepan c y in the normalization from the evolved 
IWillott et al.l (|2010[) model; we will discuss this point 
further in § 16.61 

The empirical QLF model of iHopkins et al.l (|2007bL 
hereafter HRH07) combines observations in the optical. 
X-ray, and mid-infrared bands to construct a bolomet- 
ric QLF from z = to 2; = 5. At high redshift the 
most constraining data in HRH07 comes from optical 
surveys, which have sho wn a fia ttenin g of the bright 



end slo pe at high re dshift ( S chmidt et al.il 19951 iFan et ahl 
l2001btlRichards et al. 2006 ). In the HRH07 model, the 
break luminosity increases with redshift until z ~ 2 and 
then turns over, such that it is a factor of ^ 10 lower 
luminosity at z = 5 than at z = 2, contrary to our 
finding that the break luminosity increases monotoni- 
cally to z ~ 5. At z = 5, the HRH07 model pre- 
dicts a quite low break luminosity {AI*^^q ^ —22.6) 
and a shallow bright end slope (/3 ~ —2.5). Figure [Tj] 
shows that the HRH07 model agrees with our data at 
—27 < Afi45o < —26, which is not surprising since at 
z ~ 5 their fit is mainly to the .Richards et al- (,2006i) 
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Fig. 16. — Evolutionary models for the QLF plotted against 
our binned data at 2 = 4.9. The dark and light blue dashed lines 
are the BOSS LEDE and modified LEDE forms, respectively, as in 
Fig. 1151 The LEDE model clearly disagrees with the data; the mod- 
ified form agrees with the bright end reasonabl y well. Th e black 
dotted line represents the fi t at 2 = 6 from W illott et al.l II2010I ) 
shifted to 2 = 4.9 using the IFan et al.l |j2001b ) redshift evolution; 
the shape is roughly similar to our data but it clearly understi- 
mates the no rmalization. The green dot-dashed line represents the 
QLF fit from lFontanot et al.l l(2007l '). 

data. The agreement at lower luminosities is likely a co- 
incidence due to the faint-end slope being shallow (and 
similar to their bright-end slope). However, the HRH07 
model clearly fails to match the shape at the bright end, 
which we attribute to a steeper slope and a much brighter 



Mj*45Q. This demonstrates the substantial degeneracies 
in the QLF parameters, as the HRH07 values for $* and 
Mj*45Q a,t z — 5 would not even appear on Figure [151 even 
though the model itself provides a good fit to our data at 
^^1450 ^ —27. The key diffe rence between our wo rk and 
HRH07 (and by extension, [Richards et alJ [20061 ) is the 
increased survey area at the bright end, which is needed 
to extend above the break luminosity at this redshift and 
make the break in the luminosity function more evident. 

Finally, we examine the z ~ 5 QLF from 
( |Fontanot et al.l l2007| ). This work combines bright 
quasars from SDSS DR3 with faint quasars from the 
Great Observatories Origins Deep Survey (GOODS; 
[Dickinson et al.|[2003() over the redshift range 3.5 < z < 
5.2. Figure [16] compares their pure density evolution 
(PDF) model ( #12 from their Table 3) to our data. 
Wc find that the iFontanot et alJ (|2007[ ) model overesti- 
mates our QLF measurements at all luminosities, possi- 
bly due to the comple teness correction they applied (see 
also llkeda et al.ll2012D . 

In conclusion, we find tentative evidence for a mild 
steepening of the faint-end slope, and no evidence in fa- 
vor of an evolution in the bright end slope, although both 
of these parameters are poorly constrained at most red- 
shifts. On the other hand, we do see evidence for strong 
evolution in the break luminosity, as it brightens from 
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-25.4 at z = 2.5 to M 
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-26.8 at z = 5 

(Figure [TS]) . This evolution has consequences for surveys 
where the faint limit is near the break luminosity, as sin- 
gle power law fits to such data would naturally find a 
flatter slope. The problem is evident in Figure [T^J which 
shows that a si ngle power law can d escribe the SDSS 
DR3 data from [Richards et al.l ([2006D . the fufl range of 
which is near the break luminosity. The possibility that 
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high redshift fits to the bright end of the QLF may be 
biased by a higher brea k luminosity was put forward by 
lAssef et aLI poTll) and lShen fc Kdlvl ([2012I) . Based on 
our fits at z = 5 and the greater dynamic range of our 
survey, we find this scenario to be plausible; i.e., the 
fiattening of the bright en d slop e at high redsh ift re- 
ported previously ([S chmidt et al.lll995l:lFan et al.ll2001bl: 
iRichards et"all [200& .Hopkins et al.l l2007bD imv be at- 
tributed instead to rapid evolution in the break luminos- 
ity. 

6.5. Comparison to Theoretical Predictions 

Figure [T7] compares our data to var ious theoretical 
models for the QLF ai z — h. iShenl (|2009f) provides a the- 
oretical prediction for the evolution of the QLF in a cos- 
mological framework by relating the growth of SMBHs 
to the hierarchical assembly of their host dark matter 
halos. In this model, quasar activity is triggered by ma- 
jor mergers of halos. Their fiducial model reproduces the 
observed QLF at 0.5 < z < 4.5, but underpredicts the 
observed QLF at higher redshifts. Indeed, their fi ducial 
model lies below our data at z = 5 (Figure nil). Ish^ 
(jlOOQ) also has a variant of this fiducial model that in- 
cludes a faster redshift evolution in the normalization 
of the scaling relation between peak quasar luminosity 
and host halo mass, a redshift evolution in the scatter 
of this relation and an increase in the upper host halo 
mass above which quasar triggering is cut off exponen- 
tially in order to better fit the QLF at higher redshifts 
(see Sec. 4.4 of iShenii200a . for details). This alternative 
model provides a good fit to our data at the bright end, 
b ut overpredicts th e numb er counts at the faint end. 

iConrov fc Whitd ()2012l ) introduce a model in which 
quasars are tied to galaxies in a straightforward manner 
through the Mbh - Afgai relation, so that the evolution 
of the QLF is determined by the evolution of this re- 
lation. Their model includes two free parameters that 
are allowed to vary with redshift: the quasar duty cy- 
cle and th e normalization of the A/bh ~ -Mgai relation. 
IConrov fc' Whitc (2012) constrain these parameters from 
existing QLF measurements at 0.5 < z < 4.75. The 
high redshift constra ints mainly come from the SDSS 
(|Richards et al.ll2Q06[ ) and have limited dynamic range, 
leading to significant uncertainties in the z ~ 5 predic- 
tion (represented b y the gray shaded regio n in Figure [T7l 
see also Fig. 3 in Con rov fc White! I2012[ ). Overah, the 
IConrov fc White (2012,) model provides a good fit to the 
d ata. 

iDegraf et al.l (|2010[ ) present a QLF prediction based on 
hydrodynamic simulations that include radiative cooling, 
star formation, black holes, and feedback processes. This 
work has since been updated to a larger simulation vol- 
ume, with 100 Mpc on a side and 2x1792'^ particles 
(DeGraf et al, in prep.). At each timestep, the QLF 
is calculated from the active black holes, and the final 
QLF prediction is derived by time-averaging the individ- 
ual measurements. This allows the bright end of the QLF 
to be estimated by catching the brief episodes of peak lu- 
minosity among the rare, massive black hole population. 
Figure [17] shows the prediction of this model in our red- 
shift range, including an estimate for cosmic variance de- 
rived by comparing two simulation volumes (represented 
by the extent of the shaded region). This model gener- 
ally agrees with our data near and just below the break 



10-' r 



3 10-' 
10 
10-1' 




Shen09 
CW12 

MassiveBlackll 



I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 



-21 -22 -23 -24 -25 -26 -27 -28 

^1450 

Fig. 17. — Comparison of ob served QLF to various theoretical 
predictions. The two models of IShenl l|2009l ) are shown as dashed 
brown lines; the lower line is their fiducial model, and the upper 
line is their modified m odel as described in the text. The model of 
IConrov fc Whit3 l|2012) is represented by a solid gray line, with the 
shaded region indicating the la uncertainty range. The range of 
the predicted QLF from the "MassiveBlackll" hydrodynamic sim- 
ulations is represented by the cyan shaded region. 

luminosity. It appears to be somewhat steeper than the 
data at faint luminosities, though fainter measurements 
of the QLF from deeper optical data or at other wave- 
lengths (e.g.. X-rays) are needed to better constrain the 
model. 

We remind the reader that our QLF only accounts 
for unobscured. Type I quasars. Our measurements are 
lower limits on the true density of actively accreting black 
holes, assuming some fraction are in an obscured (or even 
mildly extincted) growth phase. For example, compari- 
son to X-ray surveys indic ates that only ^ 25% of quasars 
are unobscured at z = 4 (jMasters et al.|[2"012D . Some in- 
consistency with theoretical models that do not distin- 
guish between unobscured and obscured quasars is thus 
expected. 

6.6. Spatial Density of Luminous Quasars 

A rapid decline in the comoving number density of 
quasars at high redshift was observed three decades ago 
(jOsmcr 1982). Following Fan et al. (2001a), we quantify 
this evolution in terms of the spatial density of quasars 
above a minimum luminosity within a redshift window. 
The density is derived from the 1/Va method, where for 
each quasar the volume within which it would have been 
observed within a survey is 



Az 



p{Mi4:50,z)^dz , 
az 



where p{M, z) is again the selection function for the sur- 
vey. From this equation, the total spatial density and its 
uncertainty are estimated by 
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where the sum is over all quasars more luminous than 
M. This density is related to the luminosity function in 
that 

p{< M, z)^ $(M, z)dM , (3) 
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where p(< M, z) is the space density of quasars more 
luminous than M. We choose to perform a sum over 
our data rather than an integration over the QLF as the 
latter requires extrapolation and model fitting. We cal- 
culate this quantity at 2 ^ 4, 2 ^ 5, and 2 6, com- 
bining data from SDSS, our work, and multiple surveys 
at 2 ~ 6. We choose a limit of M1450 = —26 as it cor- 
responds to the lowest luminosity quasars in the SDSS 
sample at 2 ~ 4. 

At 2 = 4.25 we use the uniform quasar sample from the 
SDSS DR7 described in i ieXTl This sample includes 311 
quasars with 4.1 < 2 < 4.4, representing a factor of four 
increase in number ov er the SDSS DR3 results given in 
[Richards et al.l (|2OO60 . We adopt the selection function 
given in Table 1 of that work rather than recalculate it 
from our simulations, as the SDSS quasar target selec- 
tion has a complicated dependence on the exten t of the 
stellar locus in col or space (.Rich ards et al] l2002l) . which 
is captured by the [Richards et a l. (2006) selection func- 
tion. At 2 = 4.9 we use our combined sample from the 
SDSS DR7 and from Stripe 82 and the selection func- 
tions presented in section 16. 1[ Finally , at 2 = 6 we use 
the sample compiled by Willott et al.l (120101) . consisting 
of quasars fr om the SDSS main (Fan et al. 2001a, 120031 
[200I I2006bn. SDSS deep (Ijiang et al.. .2008. .2009ft. and 
CFHTQS (|Willott et al.lboiOD . We derive selection func- 
tions for the three 2: 6 quasar surveys using our quasar 
simulations, extended to higher redshift. The derived se- 
lection functions from our models agree well with those 
shown in Figure 4 of Willott et al.l (|2010[ ). 

Figure [TH] shows the evolution of the space density of 
luminous quasars (M1450 < —26) from 2 = 4 to 2 = 
6. Our measurement at 2 = 5 shows a decl ine in the 
space density by a factor of 1.8 from 2 = 4.25. iFan et aTl 
([200 lb) fit an exponential decline to the space density at 
high redshifts, finding that p(Mi45o < -25.5, 2) ~ 10*"^ 
with k = —0. 47 at 2 > 3 . 6, abo ut a factor of three per 
unit redshift. iBrusa et al.l (|2009D reported a similar value 
(A; = —0.43) using X-ray-selected quasars at 2 ^ 3.0-4.5 
{\ogLx > 44.2). Our results from 2 = 4.25 to 2 = 5 
are in good agreement with these results; we measure a 
decline of 2.4 ± 0.3 per unit redshift (fc = -0.38). 

On the other hand, the evolution to 2 = 6 is more pro- 
nounced: a factor of 5.1±1.5 per unit redshift (fc = —0.7), 
or roughly twice the rate measured at lower redshift 
(Figure ITS)) . The departure of the 2 = 6 measure- 
ment from the prediction based on the slope fit to the 
4 < 2 < 5 data is significant at the 4.5(7 level. If we inte- 
grate deeper for the two high redshift bins, to a limit of 
-^1450 < —24.5, the decline is even stronger. We obtain 
p(Mi45o < -24.5, 2 = 4.9) = 4.2 ± 0.7 X 10"^ Mpc"^ and 
p(Mi45o < -24.5,2 = 6.0) = 0.5 ± 0.1 x lO'^ Mpc-^ 
corresponding to a factor of 6.5 ± 1.7 decline per unit 
redshift. 

We check the consistency of our faint-end mea- 
surements by comparing to recent work at X- 
ray wavelengths w ith the Chandra-COSMOS survey 
(jCivano et al.l 120111) , which is sensitive to moderate- 
luminosity (logL(2_io kcV) > 44.15 erg s~^) AGN at 
2 ~ 5. In order to compare the two samples, we must 
integrate our double power law mo del to Mi45n = — 24 
(adopting the Qqx relation from lYoung et al.l I2O1O0 . 
iCivano et ahl ()20lH) identified three X -ray-selected Type 
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Fig. 18. — Evolution of the number densit y of luminous quasars 
at high redshift. The data are described in i|6.6l The yellow solid 
line shows the slope derived from a log-linear fit to the 2 = 4 and 
z = 5 p oints [logp(2) oc —0.38^], a value that is similar to the 



IFan et a l. (2001b) result based on 39 bright quasars at 3.6 < 2 < 5 
from early Stripe 82 data. The 2 = 6 point is well below this 
trend, departing from it at 4.5cr. From 2 = 5 to 2 = 6 the number 
density declines as log p(2) = —0.72. The upper dark blue dashed 
line shows the extrapolated LEDE QLF model from the BOSS 
DR9, which clearly overpredicts the high redshift number densities. 
We modify the LEDE model to have a steeper evolution in logfl"* 
and shallower evolution in MJ'^gg with redshift (light blue dashed 
line, see details in text), and further impose a maximum break 
luminosity of AfJ'^gg > —27 (light blue dot-dashed line). 

1 quasars in the redshift bin 4.2 < 2 < 5.4, which 
translates to an observed space density of $(4.2 < 2 < 
5.4) = (2.5 ± 1.5) X lO"'^ Mpc"^. Integrating our best- 
fit model for the QLF over that redshift bin, we obtain 
$(4.2 < 2 < 5. 4) = 0.8 X 1(3-^ Mp c"'^. Thus our QLF 
agrees with the ICivano et al.l (|2011| ) results to within the 
~ 1(7 uncertainties, especially when considering that to 
make the comparison we must extrapolate our QLF in 
both luminosity and redshift. The X-ray data do not 
strongly constrain the evolution to 2 ^ 6 due to the small 
area of th e Ciiandra-CO SMOS survey (from zero 2^6 
objects in ICivano et al.ll201 1. the decline from 2 = 5 to 

2 = 6 is > 2 per unit redshift). 

We now relate the decline in the high-redshift quasar 
density to the evoluti onary model for t he QLF derived 
from the BOSS DR9 (|Ross et al.ll20Tl . In this LEDE 
model, log $* declines by 0.6 dex per unit redshift, while 
M* brightens by 0.68 mag per unit redshift (ci = —0.6 in 
equation [1] and C2 = —0.68 in equation [21 respectively). 
Integrating this model (equation ^ shows that it signifi- 
cantly overpredicts the high redshift number densities we 
have derived (Figure [T5)) . This result suggests that the 
steeper decline in the high-redshift number counts begins 
at 2 < 4, and is also consistent with Figure [151 which 
shows that the log predicted by the LEDE model is 
higher than the data at 2 > 5, and is not well con- 
strained at 2 ~ 4. We thus modify the LEDE model 
from BOSS by steepening the slope of the log$*(2) evo- 
lution to ci = —0.7, and softening the slope of the Mi^i^q 
evolution to C2 — —0.55. These modifications are within 
the 1(7 uncertainties of the values derived from fitting 
BOSS data; Figure [THl shows that they bring the QLF 
model in agreement with the integral constraints from 
the 2 ~ 4-5 data, but not the z = 6 point. At z = 4.9, 
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this modified evolutionary model predicts M^^^^q ~ —26.8 
and log$* = —8.0. This provides a good match to 
the values we obtained when fixing the slopes to the 
BOSS values during the MLE fit, Ml^^^ = -26.6 and 
log$*(z = 4.9) = -7.9 (Table O. 

Finally, we consider further modifying the model for 
the evolution of M* . It is unlikely that M* contin- 
ues to rise to very high luminosities; indeed, Figure [15] 
shows that the high-rcdshift fits for M* are somewhat 
below the LEDE prediction. We thus impose a maxi- 
mum break luminosity of Afj*45Q — —27, which in the 
modified LEDE model is reached at z ~ 5.4. This choice 
causes a turnover in the high-redshift number densities 
that matches the data at z 6 (Figure[T8|. This model is 
somewhat arbitrary, but qualitatively, a model in which 
$* has a more rapid downward evolution at high red- 
shift and M* brightens until z ~ 4-5 and then levels 
off (or even turns over) provides a reasonable description 
of the high redshift data. A quasar with Mj*4gQ = —27 
radiating at the Eddington luminosity corresponds to a 
~ 2 X 10^ Mq black hole; it seems reasonable to suppose 
that while the break luminosity evolves strongly at high 
redshift, it would not greatly exceed this value. 

6.7. Contribution of z ^ 5 quasars to the ionizing 
background 

The number of ionizing photons required to main- 
tain hydrogen ionization as a function of redshift can be 
found by balancing the ionizing photon emissivity with 
the density of hydrogen and the rate of recombinations. 
The recombination rate depends on the clumping factor, 
^ = {''^h) / (^h)'^ ■ iMadau et al.l ()1999f ) present an equa- 
tion for the ionizing photon density TVion that adopts a 
high value for the clumping factor (C = 30). More recent 
work suggests that the clumping factor is not so large; 

tiiksin (2005) argues that C* w 5 and iMcQuinn et al.l 
)lli) use cosmological simulations that include the ef- 
fects of Lyman- limit systems and find C « 2-3 at z = 6. 

We calculate the ionizing photon output for each 
quasar by assuming a broken power-law SE P with an in- 
dex of = —1.7 at ultraviolet wavelengths () Teller et al.l 
l2002f ). a break at llOOA, and an index of a,y = —0.5 
above the break. Rescaling the IMadau et al.l (|19990 
equation and updating to our cosmology, we estimate 
the number of photons required to maintain full ioniza- 
tion at z = 5 to be Nion = 3.4 x 10^°{C/5) Mpc"^ s-\ 
Integrating the QLF model fits to M1450 = —20 gives 
TVq ~ 6 X 10"*^ Mpc"^ s-\ or w 17% of the number 
required (this number depends only weakly on which 
model fit is used). For C = 2, quasars provide 40% 
of the photons required for hydrogen ionization, suggest- 
ing they may have played some role in maintaining ion- 
ization at z ~ 5. However, the steeper decline of lu- 
minous quasars from z ~ 5 to z ~ 6 further reduces 
the likelihood that quasars were an important source of 
ionizing photons during the reionization epoch. Indeed, 
even when as suming a steep faint-end slope (a = —1.8), 
iWillott et al.l (|2010) find that z 6 quasars generate 
< 10% of the required ionizing photon background. This 
is also in goo d agreement with c onstraints froru deep 
X-ray surveys (jBarger et al.ll2003t iFontanot et~alll2007h 
that limit the contribution from moderate luminosity 
AGN at z > 6. 



7. CONCLUSIONS 

This work builds on the le gacy of SPSS St ripe 82 for 
high-redshift quasar studies. iFan et aTl (j2004') first used 
coadded photometry from 5 epochs of Stripe 82 imag- 
ing to discover a single z ~ 6 quasar with zab = 20.5, 
fainter than the li mit adopted for the single-epoch imag- 
ing. Subsequently. iJiang et al.l ([2008. 2009) created deep, 
coadded images from the 50-60 epochs available at the 
completion of SPSS I/II to discover 11 quasars at z > 5.8 
to a limit of zab — 22.8. In this work, we have uti- 
lized the bluer SPSS bands to conduct a survey of z ^ 5 
quasars, a redshift where current constraints on the QLF 
are relatively weak. 

We define a sample of 92 candidates to a limit of 
ZAB = 22.0 using color selection criteria and with a well- 
defined selection function. From this sample, 73 objects 
have spectroscopic observations, and 71 are confirmed 
high redshift quasars. We then focus on the redshift 
range 4.7 < z < 5.1 where our completeness is highest. 
Using a sample of 52 quasars from our work on Stripe 82 
combined with 146 bright quasars from the SPSS PR7, 
we calculate the optical quasar luminosity function at 
z = 5. 

We have fit a double power law model to our observa- 
tional data at z = 5, and reach the following conclusions: 

• There is no clear evidence for evolution in the 
shape of the QLF from z 2 to z ~ 6. We 
find a steep bright-end slope (/3 ^ —3.5, Fig. [T2]) 
that roughly agrees with meas urements at both 
lower (e.g., BOSS 2.2 < z < 3.5jRoss et al.ll20lll 
and higher (z = 6, IWillott et alll2010l : IJiang et al.l 
[2009) redshifts. A relatively steep value for the 
faint end slope is favored (a ~ —1-7) in agreement 
with other high redshift results, but the faint end 
is not strongly constrained by our data. 

• We see the break in the luminosity function in 
our z ~ 5 data, finding that AI^4c^q{z — 4.9) — 
-26.8l^;^ for a = -1.7 (Fig. [HI). The strong co- 
variances between the power law slopes and the 
break luminosity lead to significant uncertainties in 
these quantities, but the best fit value for the break 
luminosity is higher than at z = 3.4, where the 
BOSS results find M^^^^iz = 3.4) = -26.0 ± 0.2. 

• The decline in the space density of luminous 
quasars at high redshift is greater than indicated 
by previous surveys. We find that while the de- 
cline in the integrated space density of quasars with 
A?i450 < —26 from z 4 to z ~ 5 is about a factor 
of three per unit redshift, in agreement with previ- 
ous results, the decline to z ~ 6 is nearly a factor of 
two greater (Fig. [T8|) . This suggests a much more 
rapid dropoff in luminous quasar activity at the 
highest redshifts currently probed by observations. 

• By comparing to a simple LEPE model for the red- 
shift evolution of BOSS quasars, we find that in 
addition to the steeper decline in number density, 
there is also an indication that the brightening of 
the break luminosity with redshift does not con- 
tinue indefinitely. A toy model in which the evolu- 
tion of <f>* is somewhat steeper than from a simple 
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extrapolation from the BOSS data, and in which 
the break luminosity evolves somewhat more slowly 
and peaks at Mi^^q = —27, provides a good match 
to the high redshift data. 

• Our model for the QLF at z ^ 5 predicts that 
quasars contribute ~ 20-50% of the ionizing pho- 
tons required to maintain hydrogen ionization at 
this redshift, with the uncertainty dominated by 
our lack of understanding of the clumping factor. 

The strong evolution in the high redshift quasar num- 
ber density we have found has implications for quasar 
surveys at even higher redshifts. Early forecasts for 
UKIDSS a nticipated r oughly 10 2 ~ 7 quasars from 
the survey (jWarren fc H ewctt 2002), m uch g reater than 
the one found to date (jMortlock et ahl [20111 ) from over 
half the survey area. Using a redshift evolution of 
log/9 oc — 0.47z ()Fan et al.ll2001bl ). extrapolation of our 
QLF predicts that there are ~ 9 quasars at6.5<z<7.5 
in the UKIDSS D RSplus release to a lim it oi Y = 19.6 
(based on Fig. 7 of lMortlock et al.ll2012D. T h e shal lower 
QLF measured at z = 6 hv iWillott et all (|2010[ ) pre- 
dicts ~ 6 quasars. Based on our results, these estimates 
should be revised downward to ^ 6 and ~ 3, respectively. 
Assuming the selection efficiency is < 50%, this brings 
the predicted number counts closer to the actual obser- 
vations. Further, the decline may continue to steepen 
at z > 6, so that the lower yield from UKIDSS may be 
due in part to an even stronger downward evolution with 
redshift. 

Significant progress has been made over the last decade 
in measuring the evolution of the high redshift quasar 
population. Nonetheless, the observations are not yet 
strongly constraining of models that make, for exam- 
ple, different predictions for the evolution of the power 
law slopes. Ongoing high redshift quasar surveys will 
improve these constraints. We are currently extending 
our work to fainter luminosities at z ~ 5 using deeper 
optical imaging data, obtaining improved measurements 
of the faint end slope and break luminosity. We will 
also use Stripe 82 data to measure the QLF at z ^ 4 
over a wide range of luminosities. Finally, Pan-STARRS 
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(jKaiser et al.ll2002l ) is obtaining shallow optical and near- 
infrared imaging, including the y-band, over an area of 
sky over that is more than twice that of the SDSS. To 
date, o ne z ^ 6 quasar has b een discovered from this 
survey (jMorganson et aLll2012| ): future work should re- 
duce the uncertainty on the bright end of the QLF at 
z ~ 6 and better constrain the strong evolution we have 
measured from existing data. 
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Fig. 19. — Spectra of z > 4.7 Stripe 82 quasars from BOSS, smoothed with a 5 pixel boxcar. In addition to the coordinates and redshift, 
each plot gives the plate-MJD-fiber designation for the BOSS observation. 
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Fic;. 20. — Spectra of Stripe 82 quasars observed with Magellan. The spectra have been rebinned to a resolution of i? f» 1200 and then 
smoothed with a 5 pixel (~ 10 A) boxcar. The first object is a repeat observation of an object with a BOSS spectrum. 
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Fig. 21. — Spectra of Stripe 82 quasars observed with MMT Red Channel, smoothed with a 7 pixel boxcar. The MMT spectra axe at low 
resolution {R ~ 400-600) and from short, 10-15m. observations, thus they have relatively lower S/N. We present the spectra with with the 
nominal flux calibrations obtained from standard star observations; however, the absolute flux calibrations have considerable uncertainties. 
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Fig. 22. — Spectra of Stripe 82 quasars observed with MMT Red Channel, continued. 
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Fig. 23. — Spectra of Stripe 82 quasars observed with MMT Red Channel, continued. 
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TABLE 4 

Stripe 82 high redshift quasar catalog 



RA (J2000) 


Dec (J2000) 


SAB 




''AB 






iAB 










JAB 


Mi450 




z 


p 


iNotes 


20: 


:39: 


:57 


.78 


+00: 


;45: 


;19, 


,4 




23, 


.37 ± 0, 


,08 


21 


.80 ± 0, 


.03 


21 


.63 ± 0, 


,08 




-24, 


,42 


4, 


,97 


0.46 


20120528 


20: 


:47: 


:23 


.18 


-00: 


;50: 


;21 


.9 




22, 


.67 ± 0, 


.04 


21 


.34 ± 


.02 


21 


.08 ± 0, 


,05 




-24, 


,79 


4, 


,89 


-1 


20110624 


20: 


:54: 


;13 


.18 


+00: 


;35: 


;03, 


.0 




22, 


.83 ± 0, 


.05 


21 


.59 ± 


.02 


21 


.46 ± 0, 


,07 




-24, 


.50 


4, 


,86 


0.74 


20120528 


21 


:10: 


;18 


.94 


+00: 


;22: 


;08, 


.8 


24.99 ± 0.26 


22, 


.79 ± 0, 


.05 


21 


.32 ± 


.02 


21 


.14 ± 0, 


,06 




-24, 


.63 


4, 


,66 


-1 


20110623 


21: 


:11: 


:58 


.01 


+00: 


;53: 


;02, 


.6 




22, 


.76 ± 0, 


.05 


21, 


.00 ± 0, 


.02 


20, 


.80 ± 0, 


,04 


21.16 ± 0.47 


-25, 


.29 


4, 


,98 


-1 


20110611 


21: 


;12: 


:25, 


.39 


-00: 


;01: 


;41, 


.3 


24.93 ± 0.23 


22, 


.42 ± 0, 


.03 


20, 


.98 ± 0, 


.01 


20, 


.81 ± 0, 


,04 




-24, 


.99 


4, 


,67 


-1 


20110623 


21: 


;14; 


:35, 


.27 


+00: 


;17: 


;33, 


,4 


25.81 ± 0.53 


22, 


.64 ± 0, 


.04 


21, 


.11 ± 0, 


.01 


20, 


.93 ± 0, 


,05 


20.32 + 0.06^= 


-24, 


.89 


4, 


,75 


0.84 


20111002 


21: 


;15; 


:47, 


.22 


+00: 


:31; 


;23, 


.8 


24.97 ± 0.27 


22, 


.40 ± 0, 


.03 


20 


.83 ± 


.01 


20 


.68 ± 0. 


.04 




-25, 


.27 


4, 


.84 


0.77 


20120528 


21: 


:18; 


:12 


.73 


-00: 


:03; 


;50, 


.7 




23, 


.72 ± 0, 


.10 


21, 


.60 ± 


.02 


21, 


.18 ± 0. 


.05 




-25, 


.64 


5. 


.46 


-1 


20120528 


21: 


;30; 


:08 


.93 


+00: 


;26; 


;09. 


.8 


26.41 ± 0.79 


22, 


.41 ± 0, 


.03 


20, 


.21 ± 0, 


.01 


20 


.48 ± 0. 


.03 




-26, 


.08 


4, 


,95 


0.62 


DR7 


21 


;43; 


:00 


.53 


-00; 


:20; 


:24, 


.6 


26.18 ± 0.68 


23, 


.58 ± 0, 


.09 


21 


.94 ± 0, 


.03 


21, 


.93 ± 0. 


.11 




-24, 


.11 


4, 


,84 


0.28 


20120527 


21: 


:50: 


:55 


.14 


-01: 


;04: 


;39, 


.8 




23, 


.15 ± 0, 


.07 


21 


.52 ± 


.03 


21 


.55 ± 0, 


,08 




-24, 


.45 


4, 


,73 


0.79 


20120527 


21 


:51: 


:51 


.83 


-00: 


;17: 


;11 


,1 


24.84 ± 0.27 


22, 


.44 ± 0, 


.04 


20 


.75 ± 


.01 


20 


.55 ± 0, 


,03 


20.13 + 0.05'= 


-25, 


.23 


4, 


,66 


-1 


20111002 


21 


;54: 


;45 


.28 


+00: 


;04: 


;21 


,7 




23, 


.57 ± 0, 


.09 


21 


.91 ± 


.03 


21 


.71 ± 0, 


,11 




-24, 


.05 


4, 


,74 


0.41 


20120528 


22: 


;00: 


:08 


.66 


+00: 


;17: 


;44, 


.9 


24.11 ± 0.10 


20, 


.60 ± 0, 


.01 


19 


.05 ± 


.00 


19 


.19 ± 0, 


,01 


19.27 + 0.07 


-27, 


.11 


4, 


,82 


0.88 


DR7 


22: 


;11: 


:41 


.01 


+00: 


;11 


;18, 


.9 




24, 


.31 ± 0, 


.17 


21 


.88 ± 


.03 


21 


.83 ± 0, 


,12 




-24, 


.84 


5, 


,23 


-1 


20110624 


22: 


;12: 


:51 


.50 


-00: 


;42: 


;30, 


,7 


24.96 ± 0.25 


21 


.64 ± 0, 


.02 


19 


.88 ± 


.00 


20 


.03 ± 0, 


,02 


20.55 + 0.25 


-26, 


.43 


4, 


,95 


0.61 


DK7 


22: 


;16: 


:44 


.01 


+00: 


;13: 


;48, 


,1 




21 


.94 ± 0, 


.02 


20 


.34 ± 


.01 


20 


.18 ± 0, 


,02 


20.25 + 0.17 


-26, 


.07 


5, 


,03 


0.40 


4200-55499-975 


22: 


;18: 


:51 


.74 


-01 


;02: 


;59, 


.0 


25.97 ± 0.86 


22, 


.37 ± 0, 


.03 


20 


.79 ± 


.01 


20 


.68 ± 0, 


,03 


20.83 ± 0.31 


-25, 


.23 


4, 


,74 


0.84 


4201-55443-370 


22: 


:19: 


:41, 


.90 


+00: 


;12: 


;56, 


,2 


25.37 ± 0.33 


22, 


.93 ± 0, 


.05 


21, 


.50 ± 0, 


.02 


21, 


.29 ± 0, 


,06 


21.05 + 0.10^= 


-24, 


.37 


4, 


,30 


-1 


20111002 


22: 


;20; 


:18, 


.48 


-01: 


;01; 


;46, 


.8 




24, 


.23 ± 0, 


.16 


21, 


.65 ± 0, 


.02 


20 


.34 ± 0. 


.02 


20.42 + 0.19 


-25, 


.83 


5. 


.62 


-1 


4201-55443-246 


22: 


:20: 


:50 


.80 


+00: 


;19; 


;59, 


.1 


23.83 ± 0.08 


21, 


.51 ± 0, 


.01 


20 


.08 ± 


.00 


20 


.04 ± 0. 


.02 


19.74 + 0.10 


-25, 


.95 


4. 


.72 


0.90 


DR7 


22: 


;22; 


:16 


.02 


-00: 


:04; 


;05, 


.6 


26.61 ± 0.98 


23, 


.81 + 0, 


.11 


21, 


.90 ± 


.03 


21, 


.83 ± 0. 


.09 


21.41 + 0.40 


-24, 


.29 


4. 


.95 


0.32 


20120528 


22: 


;25; 


:09 


.19 


-00; 


:14; 


;06, 


.8 


22.73 ± 0.03 


20, 


.41 ± 0, 


.00 


18 


.98 ± 0, 


.00 


18 


.79 + 0. 


.01 


18.86 + 0.05 


-27, 


.28 


4. 


.89 


0.72 


DR7 


22: 


:26: 


:29, 


.28 


-01: 


;09: 


;56, 


.5 




23, 


.27 ± 0, 


.07 


21 


.69 ± 


.02 


21 


.50 ± 0, 


,07 


20.76 + 0.28 


-24, 


.57 


4, 


,99 


0.55 


20120527 


22: 


:33: 


:27, 


.64 


-01 


;07: 


;04, 


,4 




23, 


.62 ± 0, 


.09 


21 


.52 ± 


.02 


21 


.49 ± 0, 


,07 




-24, 


.97 


5, 


,11 


-1 


20110624 


22: 


:36: 


:53 


.26 


+00: 


;26: 


;03, 


,0 


26.16 ± 0.66 


23, 


.63 ± 0, 


.09 


22 


.02 ± 


.03 


22 


.10 ± 0, 


,11 


22.10 + 0.76 


-23, 


.92 


4, 


,70 


-1 


20120528 


22: 


;38: 


:50 


.20 


-00: 


;27: 


;01 


,8 


26.06 ± 0.60 


23, 


.11 ± 0, 


.06 


21 


.49 ± 


.02 


21 


.36 ± 0, 


,06 


21.62 ± 0.17^= 


-25, 


.13 


5, 


,17 


-1 


20111002 


22: 


:38: 


:59 


.79 


+00: 


;56: 


;23, 


.6 


24.74 ± 0.61 


22, 


.88 ± 0, 


.05 


20 


.91 ± 


.01 


20 


.96 ± 0, 


,04 


21.43 + 0.51 


-25, 


.23 


4, 


,88 


-1 


4204-55470-710 


22: 


:39: 


:07 


.56 


+00: 


;30: 


;22, 


.6 




22, 


.97 ± 0, 


.05 


21 


.24 ± 


.02 


21 


.19 ± 0, 


,05 


21.04 + 0.34 


-25, 


.23 


5, 


,09 


0.33 


20110613 


22: 


;39: 


;25 


.51 


+00: 


;43: 


;37, 


,7 


26.10 ± 0.65 


23, 


.23 ± 0, 


.06 


21 


.88 ± 


.03 


21 


.81 ± 0, 


,08 




-24, 


.11 


4, 


,53 


-1 


20120528 


22: 


;42: 


:06 


.52 


+00: 


;25: 


;23, 


.9 


25.41 ± 0.33 


23, 


.39 ± 0, 


.07 


21 


.85 ± 


.03 


21 


.77 ± 0, 


,10 


21.73 + 0.82 


-24, 


,14 


4, 


78 


0.52 


20111002 


22: 


;45: 


:24, 


.27 


+00: 


;24; 


;14, 


,2 


26.46 ± 0.92 


23, 


.75 ± 0, 


.11 


21 


.72 ± 0, 


.02 


21 


.86 ± 0, 


,11 


21.30 ± 0.12^= 


-24, 


.88 


5, 


,16 


-1 


20111002 


23: 


:04: 


:13, 


.47 


-01; 


;11; 


;49. 


,2 


25.29 ± 0.37 


22. 


.78 ± 0. 


.04 


21, 


.24 ± 0. 


.02 


21, 


.14 + 0. 


.05 




-24. 


.75 


4. 


.75 


0.84 


20111002 


23: 


;09: 


:16, 


.76 


+00; 


;10; 


;02. 


.6 


26.31 ± 0.68 


22. 


.62 ± 0, 


.03 


20, 


.91 ± 0, 


.01 


20, 


.99 ± 0. 


,05 


21.16 + 0.27 


-25, 


.08 


4. 


.73 


0.83 


4208-55476-896 


23: 


;10: 


:00, 


.37 


-00; 


;43; 


;27. 


.0 


25.48 ± 0.33 


22. 


.62 ± 0. 


.03 


20, 


.88 ± 0, 


.01 


20, 


.87 + 0. 


.04 




-25, 


.21 


4. 


.83 


0.77 


20111002 


23: 


;12; 


:16 


.44 


+01; 


;00; 


;51, 


.5 


25.60 ± 0.44 


22, 


.71 ± 0, 


.04 


20 


.87 ± 0, 


.01 


20 


.75 ± 0. 


,03 




-25, 


.61 


5. 


,08 


-1 


4209-55478-714 


23: 


:20: 


:06, 


.60 


-00: 


;18; 


;22, 


,3 


26.08 ± 0.54 


23, 


,48 ± 0, 


,07 


21 


.58 ± 0, 


.02 


21 


.73 ± 0, 


,08 


21.25 + 0.39 


-24, 


,44 


4, 


,80 


0.77 


20111002 


23: 


:27: 


:41 


.35 


-00: 


;28: 


;03, 


,9 




23, 


,64 ± 0, 


,08 


21 


.80 ± 0, 


.02 


21 


.84 ± 0, 


,08 




-24, 


,17 


4, 


75 


0.59 


20111002 


23: 


:32: 


:23, 


.57 


+01: 


;02: 


;53, 


,5 




22, 


,83 ± 0, 


,04 


21 


.12 ± 0, 


.01 


20, 


.99 ± 0, 


,04 


20.66 + 0.23 


-24, 


,95 


4, 


,83 


-1 


4212-55447-680 


23: 


:34: 


:28, 


.47 


-00: 


;52: 


;07, 


,7 




23, 


,80 ± 0, 


,10 


21 


.99 ± 0, 


.03 


22, 


.18 ± 0, 


,11 




-23, 


,99 


4, 


77 


0.25 


20120527 


23: 


:34: 


55, 


.05 


-00; 


;10; 


;22, 


,1 




23, 


,56 ± 0, 


,07 


21 


.82 ± 0, 


.02 


21 


.33 ± 0, 


,05 


20.85 + 0.38 


-24, 


,65 


5, 


11 


-1 


20111002 


23: 


:42: 


06, 


.94 


+00; 


;36; 


;14, 


,1 




23, 


,30 ± 0, 


,06 


21 


.72 ± 0, 


.02 


21 


.74 ± 0, 


,08 


21.39 ± 0.13^= 


-24, 


,25 


4, 


74 


0.69 


20111002 


23 


:46: 


01 


.55 


-00; 


;38; 


;55, 


,2 




23, 


,64 ± 0, 


,09 


21 


.69 ± 0, 


.02 


21 


.56 ± 0, 


,07 




-24, 


,49 


4, 


,93 


0.62 


20111001 


23 


:47: 


30 


.56 


+00: 


;23: 


;06, 


,3 


25.88 ± 0.41 


22, 


,72 ± 0, 


,04 


21 


.09 ± 0, 


.01 


20, 


.75 ± 0, 


,04 




-24, 


,89 


4, 


71 


0.81 


20111002 


00: 


;05: 


;52, 


.33 


-00; 


;06; 


;55, 


.6 




25, 


.42 ± 0, 


.42 


23, 


.05 ± 0, 


.07 


20, 


.47 ± 0, 


,02 




-24, 


.42 


5, 


,86 


-1 


4216-55477-019 


00: 


;07: 


:49, 


.16 


+00; 


;41; 


;19, 


.5 


25.94 ± 0.54 


21, 


.46 ± 0, 


.01 


20, 


.03 ± 0, 


.00 


19, 


.93 ± 0. 


.01 


20.34 + 0.12 


-26, 


.10 


4. 


.83 


0.87 


DR7 


00: 


;15: 


:29, 


.86 


-00; 


;49; 


;04. 


.3 


25.57 ± 0.36 


22. 


.57 ± 0. 


.03 


21, 


.01 ± 0, 


.01 


20, 


.69 ± 0. 


.03 


20.56 ± 0.19 


-25. 


.20 


4. 


.93 


0.68 


20111001 


00: 


;23: 


:30, 


.67 


-00; 


;18; 


;36. 


,4 




23. 


.10 ± 0. 


.05 


21, 


.31 ± 0. 


.01 


20, 


.82 ± 0. 


.03 


21.50 ± 0.52 


-25. 


.10 


5. 


.06 


0.44 


4220-55447-474 


00: 


;35; 


:25, 


.29 


+00; 


;40; 


;02, 


.8 


23.73 ± 0.07 


21, 


.38 ± 0, 


.01 


19 


.66 ± 


.00 


19 


.71 ± 0. 


,01 


19.83 + 0.09 


-26, 


.42 


4, 


,76 


0.94 


DR7 


00: 


:36: 


:05, 


.64 


+01: 


;03: 


;45, 


,1 


25.17 ± 0.32 


22, 


,18 + 0, 


,02 


20 


.83 + 0, 


.01 


20, 


.77 + 0, 


,03 


20.78 + 0.19 


-25, 


,17 


4, 


72 


-1 


4221-55443-864 


00: 


:54: 


21 


.42 


-01: 


;09: 


;21, 


,6 




21, 


,28 ± 0, 


,01 


19, 


.92 ± 0, 


.00 


19, 


.54 ± 0, 


,01 


19.51 + 0.09 


-26, 


,63 


5, 


,09 


-1 


DR7 


00: 


:57: 


:03, 


.20 


+00: 


;10: 


;32, 


,3 




22, 


,98 ± 0, 


,05 


21 


.45 ± 0, 


.02 


21 


.31 ± 0, 


,07 


21.12 + 0.28 


-24, 


,61 


4, 


,84 


0.78 


20111001 


01 


:08: 


29, 


.96 


-01: 


;00: 


;15, 


,6 




23, 


,67 ± 0, 


,10 


21 


.85 ± 0, 


.03 


21 


.69 ± 0, 


,07 


21.35 + 0.52 


-24, 


,63 


5, 


11 


-1 


20111002 


01 


:10: 


:10, 


.22 


+00: 


;24: 


;19, 


,5 


25.07 ± 0.21 


22, 


,14 ± 0, 


,02 


20, 


.73 ± 0, 


.01 


20, 


.66 ± 0, 


,03 


21.37 + 0.47 


-25, 


,26 


4, 


71 


0.80 


4227-55481-528 


01 


:15: 


:44, 


.78 


+00: 


;15: 


;14, 


,9 


26.10 ± 0.58 


22, 


,97 ± 0, 


,04 


21 


.47 ± 0, 


.02 


21 


.09 ± 0, 


,05 


20.59 ± 0.06^= 


-25, 


,00 


5, 


10 


-1 


20111003 


01 


:21: 


:35, 


.57 


+00: 


;12: 


;05, 


,7 


24.51 ± 0.12 


22, 


,37 + 0, 


,03 


20, 


.98 ± 0, 


.01 


21 


.02 ± 0, 


,05 


21.24 + 0.45 


-25, 


,01 


4, 


,72 


0.83 


4228-55484-774 


01 


:28: 


:20, 


.81 


-00: 


;16: 


;36, 


,9 


25.39 ± 0.27 


22, 


,65 ± 0, 


,03 


21 


.20 ± 0, 


.01 


21 


.26 ± 0, 


,05 


21.12 + 0.48 


-24, 


,79 


4, 


,74 


0.84 


20111001 


01 


;29: 


:07, 


.44 


-00; 


;28; 


;45, 


.5 




23, 


.02 ± 0, 


.04 


21 


.54 ± 0, 


.02 


21 


.34 ± 0. 


,05 


21.31 + 0.45 


-24, 


.77 


5. 


,01 


0.58 


20111001 


01: 


;57: 


:05, 


.11 


-01; 


;12; 


;48. 


.8 


25.35 ± 0.33 


23, 


.73 ± 0. 


.08 


21, 


.92 ± 0, 


.03 


21, 


.98 ± 0. 


.09 


21.19 + 0.47 


-24. 


.08 


4. 


.79 


-1 


20111003 


01: 


;58: 


:23, 


.92 


-00; 


;36; 


;36. 


.3 




23, 


.30 ± 0. 


.05 


21, 


.62 ± 0, 


.02 


21, 


.39 ± 0. 


,05 




-24, 


.50 


4. 


.89 


0.70 


20111001 


02: 


;00: 


:01, 


.00 


+00; 


;12; 


;12. 


.9 




23. 


.15 ± 0, 


.04 


21, 


.64 ± 0, 


.02 


21, 


.40 ± 0. 


,06 


21.21 ± 0.11"= 


-24. 


.31 


4. 


.68 


-1 


20111001 


02: 


;08; 


:04, 


.31 


-01; 


;12; 


;34, 


.3 




23, 


.12 ± 0, 


.05 


21 


.46 ± 0, 


.02 


21 


.21 ± 0. 


,04 




-25, 


.28 


5. 


,23 


-1 


20111001 


02: 


:10: 


:43 


.16 


-00; 


;18; 


;18, 


,4 


22.69 ± 0.02 


20, 


,43 ± 0, 


,00 


19, 


.18 ± 0, 


.00 


19, 


.09 ± 0, 


,01 


19.33 + 0.06 


-26, 


.90 


4, 


73 


0.93 


DR7 


02: 


:10: 


:56 


.67 


+00: 


;03: 


;04, 


,7 




23, 


,97 + 0, 


,10 


21, 


.60 + 0, 


.02 


21, 


.41 ± 0, 


,07 


21.17 iO.lO'^ 


-24, 


,77 


5, 


,05 


0.47 


20111003 


02: 


ill: 


:02, 


.72 


-00: 


;09: 


;10, 


,1 


24.56 + 0.13 


21, 


,64 + 0, 


,01 


20, 


.05 + 0, 


.00 


20, 


.03 ± 0, 


,01 


20.20 + 0.12 


-26, 


,19 


4, 


,92 


0.71 


4236-55479-475 


02: 


:18: 


:21, 


.08 


-00: 


;27: 


;45, 


,7 




23, 


,72 + 0, 


,09 


21, 


.65 + 0, 


.02 


21, 


.31 ± 0, 


,04 




-24, 


,84 


5, 


11 


-1 


20111001 


02: 


:30: 


06, 


.05 


+00: 


;26: 


;25, 


,1 


25.39 + 0.25 


22, 


,94 + 0, 


,04 


21 


.69 + 


.02 


21 


.46 ± 0, 


,05 


21.65 + 0.59 


-24, 


,26 


4, 


,62 


-1 


20111002 


02: 


:36: 


:33, 


.84 


-01: 


;08: 


;39, 


,1 


25.70 ± 0.46 


23, 


,34 + 0, 


,06 


21 


.29 + 


.01 


21 


.43 ± 0, 


,05 


21.89 + 0.82 


-24, 


,97 


4, 


,97 


0.64 


20111001 


02: 


:50: 


19 


.77 


+00: 


;46: 


;50, 


,4 


22.99 ± 0.03 


21, 


,10 + 0, 


,01 


19, 


.68 + 0, 


.00 


19, 


.63 ± 0, 


,01 


19.90 + 0.11 


-26, 


,43 


4, 


,80 


0.91 


4241-55450-836 


02: 


:52: 


11 


.24 


+00: 


;44: 


;31, 


,8 


24.95 + 0.21 


22, 


,55 + 0, 


,03 


21, 


.06 + 0, 


.01 


21, 


.12 ± 0, 


,04 


21.74 + 0.72 


-24, 


,95 


4, 


,76 


0.83 


20111001 


02: 


:56: 


:17, 


.73 


+00: 


;19: 


;04, 


,2 




23, 


.56 + 0, 


.09 


21, 


.73 ± 0, 


.02 


21, 


.83 ± 0, 


,10 


21.99 ±0.13'= 


-24, 


.27 


4, 


,80 


0.66 


20111001 


02: 


;56; 


;45, 


.74 


+00; 


;02; 


;00, 


.1 




23, 


.02 ± 0, 


.05 


21, 


.58 + 0, 


.02 


21, 


.54 ± 0, 


,07 


21.15 + 0.36 


-24, 


.63 


4, 


,95 


0.67 


20111002 


03: 


;03; 


;15, 


.05 


-00; 


;03; 


;47, 


.6 


24.80 ± 0.20 


23, 


.07 + 0, 


.05 


21, 


.44 ± 0, 


.02 


21, 


.40 ± 0. 


,05 


21.12 + 0.38 


-24, 


.53 


4. 


,72 


-1 


20111003 


03: 


;10: 


:36 


.97 


-00; 


;14; 


;57, 


.0 


23.01 + 0.04 


21, 


.35 ± 0, 


.01 


19, 


.80 ± 0, 


.00 


19 


.75 ± 0. 


,01 


19.81 + 0.06 


-26, 


.24 


4. 


,72 


-1 


DR7 
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TABLE 4 — Continued 



R A f T2nnn'i 


Dcr ( 12000") 










■2^AB 








P 




03:21:55.07 


+00:29:41.4 




22.83 + 0, 


.04 


21.44 + 0.02 


21.25 + 0.05 




-24.92 


5.04 


0.52 


20111002 


03:24:59.09 


-00:57:05.0 


25.86 + 0.71 


22.53 + 0, 


.03 


20.30 + 0.01 


20.37 + 0.02 


20.01 + 0.10 


-25.76 


4.77 


0.87 


20111001 


03:27:07.86 


+00:03:44.1 




23.11 + 0, 


,06 


21.45 + 0.02 


21.45 + 0.07 


21.63 + 0.35 


-24.66 


4.88 


0.75 


20111002 


03:32:31.34 


-00:02:17.3 


26.20 + 0.71 


22.94 + 0, 


.04 


21.60 + 0.02 


21.60 + 0.07 


21.90 + 0.43 


-24.35 


4.68 


-1 


20111002 


03:38:29.31 


+00:21:56.3 


25.76 + 0.46 


21.39 + 0, 


.01 


19.84 + 0.00 


19.60 + 0.01 


19.83 + 0.08 


-26.61 


5.03 


0.35 


DR7 


03:38:30.01 


+00:18:40.1 


26.60 + 1.00 


22.75 + 0, 


,04 


21.16 + 0.01 


21.03 + 0.05 


21.48 + 0.37 


-25.09 


4.97 


0.64 


20111002 


03:46:12.30 


-00:45:12.7 




22.07 + 0, 


,03 


20.49 + 0.01 


20.46 + 0.03 


20.80 + 0.39 


-25.59 


4.80 


0.83 


20111001 


03:49:59.40 


+00:34:03.5 




23.54 + 0, 


,11 


21.38 + 0.02 


21.27 + 0.06 


21.28 + 0.30 


-25.31 


5.20 


-1 


20111002 


03:50:43.92 


+00:17:17.7 


25.27 + 0.47 


22.56 + 0, 


,05 


21.28 + 0.02 


21.43 + 0.10 




-24.72 


4.75 


0.84 


20111003 



Note. — Magnitudes are on the AB system lOke fc Gunnin"983l ) and have been corrected for extinction. Blank entries indicate objects 
for which Sextractor did not report a magnitude in that band (due to negative fluxes); objects that are not formally detected but have 
measurements from Sextractor are included with their uncertainties. Ellipses indicate lack of coverage. 

^ The selection probability p{M,z) of spectroscopically confirming this object, including all sources of incompleteness (see i|6. 1.21 1. A 
value of —1 indicates an object not included in the uniform sample. 

^ Spectroscopic identifications from MMT/Magellan show the UT date of the observation; those from BOSS give the plate-MJD-fiber of 
the observation. 

J magnitude from SWIRC. 
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TABLE 5 

Stripe 82 high redshift quasar targets 



RA (J2000) Dec (J2000) ^ab ^ab «ab -zab -'ab Notes 



20; 


:37: 


:53, 


,64 


+00: 


:05: 


:49. 


,4 












22, 


.81 + 0, 


,05 


20, 


.84 + 0.01 


20. 


.33 + 0. 


.03 






20: 


:49: 


:47, 


.51 


+00: 


:08: 


:24. 


1 


25, 


.82 


± 


0, 


.59 


23, 


.38 + 0. 


,09 


21, 


.91 ± 0.03 


21, 


.64 + 0, 


.11 




star 


20: 


:59: 


:54, 


.11 


-00: 


:54: 


:27, 


.7 


25, 


.90 


± 


0, 


.76 


23, 


.11 + 0. 


,06 


21, 


.54 + 0.02 


21, 


.24 + 0, 


.05 


21.22 + 0.15=' 




21: 


:00: 


:41, 


.31 


-00: 


:52: 


:03, 


.4 


25, 


.58 


± 


0, 


.58 


23, 


.59 + 0. 


10 


21, 


.79 ± 0.03 


21, 


.37 + 0, 


.06 


20.94 + 0.10=^ 




21: 


:05: 


:17, 


.96 


-00: 


:49: 


:20, 


.2 












23, 


.23 + 0. 


,07 


21, 


.49 ± 0.02 


21, 


.18 + 0, 


.06 


20.86 + 0.31 




21: 


:47: 


:30, 


.25 


+00: 


:12: 


:37. 


,5 


25, 


.32 


± 


0, 


.38 


23, 


.65 + 0. 


11 


21, 


.92 ± 0.03 


21, 


.50 + 0. 


.09 


20.82 + 0.38 




21: 


:52: 


:56, 


.25 


-01: 


:06: 


:13, 


.8 


25, 


.01 


± 


0, 


.46 


23, 


.00 + 0. 


,07 


21, 


.56 ± 0.03 


21, 


.32 + 0. 


.07 


20.66 + 0.08=^ 




00: 


:03: 


:32, 


.25 


+00: 


:37: 


:49. 


,9 


26, 


.30 


± 


0, 


.72 


23, 


.18 + 0. 


,05 


21, 


.61 ± 0.02 


21, 


.37 + 0, 


.05 


21.30 + 0.32 




00: 


:05: 


:31, 


.71 


-00: 


:34: 


:43, 


.4 












23, 


.85 + 0. 


10 


21, 


.91 ± 0.03 


21, 


.75 ± 0, 


.07 


21.77 + 0.64 




uu. 


■1 

, lo. 


■99 


1 n 


-rUU. 


.or:. 
, ZO. 


9*^ 


1 

, 1 


ZD. 


90 
,ZU 


_|_ 
□I 


n 
u. 




9*^ 

Zo. 


fil -1- n 
.ui m u. 




91 

Zl , 


7Q _|_ n 09 

. / y III u.uz 


91 
Zl, 


.Do m u. 


.Uo 






00: 


:23: 


:53, 


.63 


-00: 


:43: 


:56, 


.9 


26, 


.00 


± 


0, 


.52 


23, 


.64 + 0. 


,08 


21, 


.51 ± 0.02 


21, 


.00 ± 0. 


.04 


20.75 + 0.07^ 




00: 


:30: 


:30, 


.36 


-00: 


:04: 


:17, 


.7 


25, 


.78 


± 


0, 


.42 


23, 


.41 + 0. 


,06 


21, 


.78 ± 0.02 


21, 


.67 + 0. 


.06 


21.69 + 0.76 




00: 


:41: 


:37, 


.17 


-00: 


:30: 


:39, 


.4 












22, 


.94 + 0. 


,04 


20, 


.92 ± 0.01 


20, 


.45 ± 0. 


.02 


20.00 + 0.15 




01: 


:26: 


:55, 


.05 


-00: 


:00: 


:12, 


.7 












24, 


.37 + 0. 


16 


21, 


.87 + 0.02 


21, 


.78 ± 0, 


.10 


21.09 + 0.40 


star 


01: 


:40: 


:14, 


.61 


-00: 


:12: 


:59, 


.9 


25, 


.42 


± 


0, 


.30 


23, 


.30 + 0. 


,05 


22, 


.00 + 0.02 


21, 


.81 + 0. 


.07 


21.47 + 0.60 




02: 


:38: 


:09, 


.21 


-01: 


:09: 


:27, 


.5 


28, 


.23 


± 


4, 


.72 


23, 


.44 + 0. 


,06 


21, 


.94 + 0.03 


21, 


.83 + 0, 


.08 






02: 


:52: 


:29, 


.91 


-00: 


:48: 


:13, 


.1 












23, 


.74 + 0. 


,09 


21, 


.97 + 0.03 


21, 


.77 + 0. 


.08 


21.14 + 0.41 




03: 


:09: 


:37, 


.66 


+00: 


:46: 


16. 


,8 


26, 


.77 


± 


1, 


.16 


23, 


.32 + 0. 


,06 


21, 


.91 ± 0.03 


21, 


.74 ± 0. 


.07 


21.38 + 0.37 




03: 


:40: 


:50, 


.33 


+00: 


:48: 


12. 


,2 


25, 


.99 


± 


0, 


.61 


23, 


.23 + 0. 


,07 


21, 


.88 ± 0.03 


21, 


.73 ± 0, 


.07 


21.06 + 0.23 




03: 


:54: 


:23, 


.45 


+00: 


32: 


31. 


,6 


24, 


.00 


± 


0, 


.44 


22, 


.92 + 0. 


15 


21, 


.28 ± 0.04 


21, 


.05 ± 0, 


.08 


20.33 + 0.18 




03: 


:56: 


:55, 


.52 


-00: 


:34: 


:47, 


.9 


24, 


.82 


± 


0, 


.65 


22, 


.48 + 0. 


.08 


21, 


.21 ± 0.03 


21, 


,09 + 0, 


.07 


22.01 ± 1.16 





Note. — Targets selected by our color criteria that either lack spectroscopic observations or are confirmed non-quasars. Magnitudes 
are given in the same format as Table |4] 



J magnitude from SWIRC. 



